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ABSTRACT 


A  model  has  been  developed  to  evaluate  the  influence  of  polarization  filtering  on 
the  detection  and  identification  range  of  a  thermal  sensor.  The  scenarios  evaluated  were 
based  on  environmental  parameters  and  ship  temperatures  recorded  during  the  EOPACE 
measurement  series  in  San  Diego  bay  in  1996.  These  scenarios  represent  a  FLIR  sensor 
on  a  platform  in  level  flight  at  100  m  or  1000  m  approaching  a  ship  target  represented  as 
a  gray  body  at  the  recorded  ship  hot-spot  temperature.  The  polarized  version  of  the 
SEARAD  sea  radiance  code  was  used  to  provide  sea  background  radiance  and 
propagation  characteristics  for  both  ship  target  and  background.  Apparent  Temperature 
Difference  was  calculated  versus  range  for  horizontally  polarized  and  unpolarized 
imaging.  Maximum  range  was  determined  for  both  cases  by  comparison  to  a  generic 
Minimum  Resolvable  Temperature  function  representing  a  typical  LWIR  Common 
Module  FLIR.  Preliminary  results  for  the  polarized  case  predict  greater  apparent 
temperature  difference  at  ranges  to  around  10  kilometers.  Unresolved  apparent  anomalies 
in  the  computed  results  suggest  that  target  temperatures  are  under-estimated.  Empirical 
correction  of  the  zero  range  temperature  difference  suggests  polarized  identification 
ranges  of  the  order  of  25  to  30  km.  Improvements  to  the  modeling  are  proposed. 
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I.  INTRODUCTION 

Infrared  systems  have  numerous  applications  in  today’s  military  environment. 
Their  primary  use  is  for  target  detection,  recognition  and  identification  by  discriminating 
between  target  and  background  temperature.  There  is  an  intensive  scientific  and 
technological  effort  to  improve  the  performance  of  such  systems  due  to  the  high  potential 
of  the  military  applications  using  this  technology. 

The  operating  environment  of  the  infrared  devices  is  a  noisy  one  and  strongly 
influenced  by  the  weather  conditions.  In  order  to  improve  the  target-to-background 
contrast,  which  eventually  will  improve  the  detection  capabilities  of  a  system,  various 
techniques  have  been  tried.  One  such  promising  technique  is  the  use  of  polarization  filters 
in  the  sensor.  The  application  of  this  method  requires  an  environment  in  which  the 
electromagnetic  energy  coming  from  the  source  or  background  is  polarized.  The  sea 
background  comprises  such  an  environment. 

It  has  been  found  that  the  infrared  electric  vector  of  sea  radiance  is  preferentially 
horizontally  polarized  in  the  sun  glint  corridor  within  a  few  degrees  of  the  horizon  and 
vertically  polarized  from  the  adjacent  sea  background  at  similar  grazing  angles.  Previous 
experiments  and  analysis  done  at  the  Naval  Postgraduate  School  [Ref. 6, 13, 16,22]  have 
shown  that  there  is  a  significant  contrast  improvement  between  a  man-made  target,  such 
as  a  ship  which  gives  no  apparent  polarized  feature  properties  to  the  reflected  or  emitted 
energy,  and  the  sea.  Therefore  horizontally  polarized  fdters  can  improve  the  target-to- 
background  contrast  while  viewing  target  scenes  outside  the  sun  glint.  It  must  be  taken 
into  consideration  though,  that  the  magnitude  of  degree  of  polarization  was  found  to  be 
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much  greater  in  the  3  to  5  /Jm  band  in  the  Sun  glint  regions  than  the  8  to  12  jjm  band 
which  dominates  emission  and  is  the  band  that  will  be  used  for  this  work. 

In  this  thesis  data  from  an  experimental  data  set  that  was  obtained  in  the  Electro- 
Optic  Propagation  Assessment  in  the  Coastal  Environment  (EOPACE)  multinational 
measurements  series  in  San  Diego  Bay  in  1996,  is  used  to  retrieve  some  data  for  the 
construction  of  the  scenarios  that  were  used  for  the  analysis.  EOPACE  is  a  program 
sponsored  by  ONR  and  organized  by  SPAWAR  SYSTEMS  CENTER,  San  Diego,  Ocean 
and  Atmospheric  Sciences  Division,  Code  543.  The  purpose  of  the  EOPACE 
measurements  is  to  quantify  infrared  propagation  characteristics  for  near  ocean  surface 
transmission  and  analyze  Electro-Optic  systems  performance  in  the  coastal  environment 
conditions.  The  data  that  were  used  in  this  thesis  were  obtained  by  the  NPS  Boundary 
Layer  Meteorology  Group  on  the  Research  Vessel  (R/V)  POINT  SUR  and  NPS  NACIT 
(Infrared  Technology  Group)  at  Building  15  of  the  Naval  Command  Control  and  Ocean 
Surveillance  Center  (NRAD)  on  Point  Loma,  California,  on  9  April  1996.  The  collected 
data  comprise  5,508  images  of  the  R/V  POINT  SUR  in  sets  of  various  ranges,  aspect 
angles,  meteorological  conditions  and  ship  skin  temperatures.  The  data  that  were 
collected  in  the  experiment  are  image  sets  of  horizontally  polarized,  vertically  polarized 
and  unpolarized  images  of  the  R/V  POINT  SUR  taken  with  an  AGA  780  camera  using 
internal  polarizers  (Aluminum  grid  on  a  KRS-5  substrate). 

This  thesis  addresses  the  prediction  of  detection,  recognition  and  identification 
ranges  with  and  without  polarization  for  a  target  with  the  characteristics  of  R/V  POINT 
SUR  by  taking  under  consideration  the  specific  environmental  conditions  when  the 
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measurements  were  conducted.  The  project  addresses  two  different  scenarios  in  which  a 
generic  Forward  Looking  Infrared  (FLIR)  platform  sensor  approaches  the  target  (R/V 
POINT  SUR)  at  different  heights.  For  the  calculation  of  the  apparent  target  to 
background  temperature  difference  (ATapp)  at  given  range  the  SEARAD  atmospheric 

propagation  code  is  used,  which  is  a  modified  form  of  the  standard  MODTRAN  code 
taking  into  consideration  the  effects  of  polarization.  SEARAD  is  a  DOS  compatible 
program  developed  by  C.R.  Zeisse  at  NCCOSC-NRaD  (Naval  Command,  Control  and 
OCEAN  Surveillance  Center-Naval  Research  and  Development),  for  application  to  naval 
problems  [Ref. 28].  The  sensor  performance  is  described  by  the  parameter  Minimum 
Detectable  Temperature  Difference  (MDTD)  or  Minimum  Resolvable  Temperature 
Difference  (MRTD)  as  appropriate.  By  comparing  the  MRTD  with  ATapp  the  maximum 

ranges  for  detection,  recognition  and  identification  will  be  estimated  by  using  specific 
(standard)  criteria  to  describe  the  actual  target. 
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II.  INFRARED  TECHNOLOGY  FUNDAMENTALS 


In  order  to  follow  a  progressive  approach  to  the  objectives  of  this  thesis  a  brief 
presentation  of  the  basic  laws  and  theories  of  Infrared  technology  is  necessary.  Only  the 
most  important  points  necessary  to  comprehend  the  basic  physical  phenomena  associated 
with  this  work  will  be  addressed. 

A.  THE  ELECTROMAGNETIC  SPECTRUM 

In  general,  waves  are  means  of  transporting  energy  or  information.  Typical 
examples  of  Electromagnetic  Waves  are  light  rays,  radio  waves  and  radar  beams.  Their 
common  characteristic  is  the  spatial  and  temporal  variations  of  the  electric  and  the 
magnetic  fields.  They  relate  with  each  other  and  their  relationship  is  presented  by 
Maxwell’s  equations: 


Faraday’s  Law: 


VxE  =  -— 
dt 


(2.1) 


Ampere’s  Law:  V  x  H  =  J  + 


dV 

dt 


(2.2) 


Gauss’s  Law:  V  •  D  =  p  v  (2.3) 

Nonexistence  of  isolated  magnetic  charge:  V  •  B  =  0  (2.4) 


where: 

E=  Electric  field  intensity  vector  (Volts/meter) 
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B=  Magnetic  flux  density  (Tesla) 

T=  Time  (seconds) 

H=  Magnetic  field  intensity  vector  (Amperes/meter) 

J=  Current  density  (Amperes/meter2) 

D=  Electric  field  density  vector  (Coulomb/  meter2) 
pv  =  Volume  charge  density  (Coulomb/  meter3) 

Since  all  of  the  above  mentioned  waves  as  well  as  infrared  waves  comply  with  the 
same  laws,  they  are  of  the  same  nature.  However  quantization  effects  which  may  be 
ignored  at  long  wavelengths  are  of  major  significance  in  the  visible  and  infrared. 

What  is  known  as  optical  radiation  covers  the  electromagnetic  spectrum  from  the 
Ultraviolet  (UV)  range  throughout  the  visible  and  up  to  the  infrared  (IR)  portion  of  the 
spectrum,  as  shown  in  Figure  1  [Ref.l]. 

The  ultraviolet  portion  ranges  from  about  0. 1  up  to  0.40  pm .  The  visible  range  of 
the  electromagnetic  spectrum  to  the  human  eye,  commonly  known  as  visible  light,  is 
from  0.40 pm  up  to  0.77 pm.  The  infrared  portion  is  divided  into  three  operational 
regions;  the  near  infrared  or  short  wavelength  infrared  (SWIR)  region  ranging  from 
0.77  pm  up  to  3  pm ;  the  middle  wavelength  infrared  (MWIR)  ranging  from  3  pm  up  to 
8  pm ;  the  long  wavelength  infrared  (LWIR)  from  8  pm  up  to  14  pm;  and  finally  the  far 
and  extreme  infrared  region  ranging  14 pm  up  to  1000 pm.  The  region  from  5 pm  to 
8  pm  is  also  a  part  of  the  infrared  spectrum,  but  is  not  used  due  to  the  high  absorption  of 
the  atmosphere  in  this  area. 
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Figure  1  -  The  Electromagnetic  Spectrum  [Ref.l]. 
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Imaging  devices  working  in  the  near  infrared  region  sense  ambient  radiation  that 
is  reflected  from  the  targets  and  the  background  scene.  These  devices  are  known  as  night 
vision  devices  and  can  detect  infrared  with  wavelength  up  to  about  1.5  flm .  Infrared 
radiation  with  wavelengths  greater  than  3  [im  is  known  as  thermal  radiation,  because  the 
energy  emitted  by  the  target  due  to  its  temperature  provides  the  largest  part  of  the 
detectable  signature.  Thermal  imaging  systems  basically  take  advantage  of  this  thermal 
energy  that  every  object  emits,  and  work  in  the  regions  of  MWIR  and  LWIR.  Such 
devices  are  the  Forward  Looking  Infrared  (FLIR)  and  the  Infrared  Search  and  Track 
(IRST). 


B.  RADIANT  ENERGY  TRANSFER  UNITS 

The  basic  radiometric  quantities  and  units  of  the  infrared  radiation  are  presented 
in  Table  1  [Ref.  2:p.l4]  in  both  their  energy  and  photon  forms.  To  distinguish  between 
the  two  a  different  subscript  is  used  (e  for  energy  and  q  for  quanta).  The  most  common 
detector  type  that  is  used  is  the  quantum  detector,  whose  output  signal  is  directly 
proportional  to  the  integral  of  the  incident  photon  flux  over  a  specific  time  period.  This  is 
the  main  factor  that  makes  desirable  the  use  of  the  photon  forms.  Historically  the  energy 
forms  were  and  still  are  being  used  for  detector  performance  measurements.  Therefore 
both  approaches  are  still  useful. 
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Table  1  -  Definition  of  Radiometric  Quantities  and  units. 


Quantity 

Definition 

Photon  Symbol 

Radiant  Energy  Q 

Total  Radiated 
Energy 

Qe(J) 

Qq  (Photons) 

Radiant 

Flux  <t> 

dQ 

dt 

<I>e  (W) 

3>9  (Photons  sec'1) 

Radiant  Flux  or 
Exitance  M 

d& 

dA 

Mq  (Photons  sec'1  m' 

2) 

Irradiance  I 

d& 

3A 

Eq  (Photons  sec'1  m'2) 

Radiance  L 

a2o 

(3A  *  cos  tfjdQ 

mm 

Radiant  Intensity  I 

3<D 

3Q 

Iq  (Photons  sec-1  sr1) 

C.  BASIC  LAWS  OF  THERMAL  RADIATION 

1.  Black  Body  Radiation 

The  definition  of  some  important  parameters  is  necessary  before  addressing  the 
basic  laws  of  thermal  radiation.  Emissivity  of  a  surface  is  defined  as  the  "ratio  of  emitted 
radiant  power  from  a  surface  to  that  emitted  from  a  black  body  (perfect  emitter)  at  the 
same  temperature"  [Ref.3].  Absorptance  is  defined  "as  the  fraction  of  incident  power 
absorbed"  [Ref.3].  The  spectral  radiant  emittance  is  the  radiant  flux  leaving  an 
infinitesimal  area  of  surface  divided  by  that  area  per  unit  wavelength  [Ref.3]. 

A  black  body  can  be  defined  as  a  perfect  absorber  (absorptivity  a- 1).  Being  a 
perfect  absorber,  it  must  also  be  a  perfect  emitter  (emissivity  £  =  1).  According  to  those 
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two  properties  and  Kirchoff’s  Law  the  Black  Body  has  a  reflectance  ( p )  and 
transmissivity  ( T )  of  zero.  A  black  body  is  totally  diffuse,  radiates  at  all  wavelengths, 
and  at  any  temperature  its  spectral  radiant  emittance  is  the  maximum  possible  for  any  real 
thermal  source  with  the  same  temperature.  In  Figure  2  the  spectral  radiant  emittance  of  a 
blackbody  at  various  temperatures  is  presented. 


Figure  2  -  Spectral  Radiant  Emittance  of  a  Blackbody  at  various  Temperatures  [Ref.3:p.3-3]. 
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2. 


Total  Power  Law 


When  there  is  incident  radiation  on  a  body,  some  of  it  is  absorbed,  some  of  it 
reflected  and  some  of  it  transmitted.  Thus  the  ratios  of  each  of  those  three  to  the  incident 
power  must  add  to  be  equal  to  one: 

a+  p+r  =  1  (2.5) 


where: 

a=  Absorptivity  (the  fraction  of  incident  power  absorbed) 
p  =  Reflectivity  (the  fraction  of  incident  power  reflected) 
r  =  Transmissivity  (the  fraction  of  incident  power  transmitted) 

3.  Kirchoffs  Law 

A  non-Black  body  in  a  thermal  equilibrium  receives  but  also  emits  power.  Since 
the  system  is  in  thermal  equilibrium,  the  absorbed  and  emitted  powers  must  be  equal.  It 
follows  by  extension  that  at  any  wavelength  the  emissivity,  the  fractional  efficiency  of 
radiation  relative  to  a  black  body,  must  equal  the  absorptivity. 

e{X)  =  a{X)  (2.6) 

This  expression  is  a  statement  of  Kirchoffs  Law  together  with  the  principle  of 
“detailed  balance”. 
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4. 


Planck’s  Radiation  Law 


In  the  late  1890’s  the  German  theoretical  physicist  Max  Planck  found  a  single 
formula  that  could  fit  the  measurements  of  MA(L) ,  known  as  the  spectral  exitance,  over 
the  entire  wavelength  range.  This  formula  became  known  as  Planck’s  Radiation  Law: 


InShltf 

=  — - 

(em  -1) 

where: 

Ma  (T)  =Spectral  Exitance  (Watts/cm2/  fJm ) 

A  =  Wavelength  ( /i  m) 
k  =  Boltzmann’s  Constant  (k=1.38xl0'23  J/K) 
c=  Light  velocity  in  Vacuum  (c=2.998xl08  m/sec) 
h  =  Plank’s  Constant  (h=6.625xl0'34  Js) 

T  =  Source  Temperature  (K) 


(2.7) 


Planck  came  to  this  result  by  making  two  important  modifications  to  the  classical 

theory: 

1.  The  emission  of  absorption  of  energy  out  of  a  surface  is  by  oscillators  (of 
electromagnetic  origin)  which  can  only  have  certain  discrete  energies 
determined  by  En  =  nhf  ,  where  n  is  an  integer,  /  is  the  frequency  in  Hertz 
and  h  is  called  Plank’s  constant  equal  to  h=6.625xl0'34  Js. 
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2.  The  oscillators  can  absorb  or  emit  energy  in  discrete  multiples  of  the 
fundamental  quantum  of  energy  given  by  A E  =  hf  . 

Planck’s  Law  represents  with  great  accuracy  the  radiation  from  a  surface,  and  with 
a  table  of  emissivities  allows  the  calculation  of  the  emittance  from  a  real  surface. 

Emissivity  is  defined  as  the  ratio  of  the  radiant  exitance  or  radiance  of  a  given 
body  to  that  of  a  black  body.  The  basic  definition  is  in  terms  of  a  narrow  spectral  interval. 
Thus  the  spectral  emissivity  of  a  body  is  defined  as: 


(2.8) 


A  material  which  has  an  emissivity  that  is  independent  of  wavelength  is  often 
called  a  gray  body.  Bodies  which  have  an  emissivity  which  varies  with  wavelength  are 
often  called  spectral  or  colored  bodies,  or  selective  emitters. 


5.  Wien’s  Displacement  and  Stefan-Boltzmann  Laws 

Blackbody  radiation  is  theoretically  interesting  because  of  its  universal  character. 
The  radiation  properties  of  the  blackbody  are  independent  of  the  particular  material  of 
which  the  container  is  made. 
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The  spectral  exitance  Mx(X,T)  is  the  total  power  radiated  per  unit  area  per  unit 
wavelength  at  a  specific  temperature.  From  the  curves  in  Figure  3  two  important 
observations  should  be  noted: 

1.  There  is  a  shift  in  the  maximum  towards  the  shorter  wavelengths  as  the 
temperature  is  increased. 

2.  The  total  radiated  power  increases  as  the  temperature  increases. 

The  first  of  the  above  mentioned  observations  is  commonly  referred  to  as  Wien’s 
Displacement  law  and  is  expressed  by  the  formula: 

X^J  =2.898x1 0‘3  (m  K)  (2.9) 

Where  Xmax  is  the  wavelength  of  the  peak  of  the  spectral  distribution  at  a  given 
temperature. 

For  the  second  observation  if  we  integrate  the  quantity  Mx{X,T)  over  all 

wavelengths,  the  power  per  unit  area  at  T  is  obtained.  Josef  Stefan  found  empirically  in 
1879  and  Boltzmann  demonstrated  theoretically  several  years  later,  that  M(T)  is  related 
to  temperature  by: 


M(T)  =  soT4  (2.10) 

This  is  known  as  the  Stefan-Boltzmann  law,  with  the  constant  <7  experimentally 
measured  to  be  5.6705xl0'8  (W/m2K4). 
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Figure  3  -  Spectral  distribution  of  the  radiation  emitted  from  a  blackbody  for  different 

blackbody  temperatures  [Ref. 2]. 


D.  ATMOSPHERIC  PROPAGATION 

The  atmosphere  is  an  important  consideration  regarding  the  performance  of  all 
infrared  systems.  An  infrared  system  can  be  described  as  having  three  basic  components: 
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the  source  of  the  radiation,  the  propagation  medium  and  finally  the  detector.  Due  to  the 
dramatic  improvement  in  the  performance  of  detection  systems  today,  the  limiting  factor 
of  the  overall  system  performance  is  often  the  propagation  medium. 

Absorption,  scattering  and  turbulence  are  the  dominant  mechanisms  of  loss. 


1.  Absorption 

The  atmospheric  window  regions  are  defined  by  molecular  absorption  (primarily 
water  vapor  and  carbon  dioxide).  The  energy  structure  of  a  molecule  is  composed  of 
discrete  energy  levels.  At  optical  frequencies,  the  quantum  nature  of  molecules  is  very 
important.  Figure  4  illustrates  the  concept  for  a  two  level  system.  A  photon  with  a 
frequency  f0,  such  that  Ei-E2=hf0  (h  is  Planck’s  constant)  is  absorbed  by  the  system.  Other 
frequencies  do  not  interact.  Energy  in  this  case  is  absorbed  by  the  molecule  (from  the 
photon  field)  at  the  frequency  fo.  Thus  an  absorption  spectral  line  is  observed  in  a  very 
narrow  linewidth  around  fo  because  of  the  absence  of  energy  relative  to  the  incident 
photon  field. 


fo 


► 


fi 


E2 

E)-E2=hfo 


fi 


► 


Figure  4  -  Photon  absorption  in  a  two  level  system. 
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Figure  5  -  Atmospheric  transmittance  for  marine  aerosol  model  using  the  1976  US  Standard 
atmospheric  model  for  a  slant  range  of  2  km  from  sea  level  to  100m  height  calculated  with  the 
SEAR  AD  code. 


The  plot  in  Figure  5  shows  the  infrared  "windows"  in  which  systems  have  to 
operate.  The  extinction  that  is  shown  is  caused  by  four  separate  mechanisms  that  must  be 
considered  separately.  The  total  extinction  coefficient  is  the  sum  of  the  coefficients  for 
total  absorption  and  total  non-forward  scattering: 


/*  =  /**  +Ms 


(2.11) 
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where: 


jua  =  total  absorption  coefficient 

/js  =  total  non-forward  scattering  coefficient. 

The  total  absorption  coefficient  is: 

H.=K+k.  (2.12) 


where: 

km  =  molecular  absorption  coefficient 

ka  =  aerosol  absorption  coefficient. 

As  for  the  molecular  extinction,  the  relative  values  of  these  coefficients  depend 
strongly  on  the  density  and  molecular  composition  of  the  atmosphere  and  the 
composition,  number  density  and  size  distribution  of  the  aerosols. 

2.  Scattering 

Scattering  is  the  process  by  which  the  energy  in  an  electromagnetic  wave  is 
intercepted  and  reradiated  into  An  steradian  solid  angle.  It  results  from  the  interaction  of 
the  wave  field  with  the  electron  oscillators  in  the  scattering  medium  [Ref.  3  :p.  12-22].  If 
the  energy  of  an  electromagnetic  wave  does  not  match  the  energy  of  one  of  the  atom’s 
energy  state  transitions,  the  atom  will  scatter  the  light.  So  in  contrast  to  the  excitation 
process  (absorption)  non-resonant  scattering  occurs  with  incoming  radiant  energy  of 
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other  frequencies  (other  than  resonant  frequencies).  So  the  atoms  (or  the  molecules)  of 
the  scattering  medium  are  excited  by  the  incoming  wave  field,  behaving  as  forced 
harmonic  oscillators  which  reradiate  at  the  frequency  of  the  incident  wave.  There  is  a 
phase  difference  in  the  secondary  wave  from  the  primary,  however  (phase  lag), 
dependent  on  the  difference  between  the  wave  frequency  and  the  oscillator  resonance 
frequency.  Any  inhomogeneity  in  refractive  index  can  cause  scattering  (in  the  atmosphere 
the  component  gas  molecules,  aerosol  particles  of  various  sizes,  fog,  rain  and  hail  drops 
are  all  effective  scatterers).  The  spatial  distribution  of  the  scattered  radiation  is  strongly 
dependent  on  the  relative  magnitudes  of  the  particle  size  and  wavelength.  For  particles 
very  small  compared  to  wavelength,  the  scattering  is  approximately  isotropic.  As  the  ratio 
of  size  to  wavelength  increases,  the  scattering  is  basically  concentrated  into  the  forward 
hemisphere.  For  very  large  scattering  objects  forward  scattering  dominates  and  secondary 
lobes  develop  at  other  angles  in  the  radiation  pattern. 

The  scattering  coefficient  that  is  part  of  the  total  extinction  coefficient 
{l  =  jua  +  Jls,  can  be  divided  into  components  due  to  the  molecules  of  the  air  and  the 
aerosol  particles  suspended  in  it: 


U  =  (7  +  <7 

r*s  m  a 


(2.13) 


where: 

<7m  =  molecular  scattering  coefficient 
<Ja  =  aerosol  scattering  coefficient. 
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The  approximate  order  of  importance  in  each  of  the  important  atmospheric 


transmission  windows  is  shown  in  Table  2  [Ref.  3:p.l2-9]. 


Table  2  -  Order  of  Significance  of  Attenuation  Coefficients  [Ref.  3:p.l2-9]. 


Atmospheric  Window 

Wavelength 

Attenuation  Coefficients  in 
order  of  importance 

Visible 

0.39  (Jm  to  0.77  /Jm 

G a  »  5  ka 

Near  Infrared  (SWIR) 

0.77  /lm  to  1.2  /Jm 

Mid  Infrared  (MWIR) 

3  /Jm  to  5  ftm 

® a  ’  ’  & m 

Far  Infrared  (LWIR) 

8  /Jm  to  12  fjm 

km  >  K  ’  °a 

a)  Rayleigh  Scattering 

In  the  small  particle  region  the  scattering  can  be  described  accurately  by 
the  relatively  simple  theory  developed  by  Rayleigh  for  small  particles.  His  theory 
predicts  that  in  this  case  scattering  is  proportional  to  the  second  power  of  the  volume  of 
the  particle  and  inversely  proportional  to  the  fourth  power  of  the  wavelength.  This 
describes  the  scattering  by  atmospheric  molecules. 

The  attenuation  of  a  light  beam  due  to  Rayleigh  scattering  is  expressed  by: 

1  =  h  exp (~<ymR)  (2.14) 


where: 

Io=  irradiance  of  the  original  radiation 
1=  irradiance  of  the  remaining  radiation 
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£ Tm  =  Rayleigh  scattering  coefficient 


R=  path  length. 


b)  Mie  Scattering 

In  1908  Gustav  Mie  published  a  rigorous  solution  of  the  scattering 
problem  for  homogeneous  spherical  particles  of  any  size.  Although  complicated,  his 
solution  has  great  practical  value,  particularly  when  applied  to  colloidal  and  metallic 
suspensions,  fog  and  clouds  to  mention  only  a  few  applications[Ref.4:p.294].  Mie 
scattering  includes  Rayleigh  scattering  in  the  limit  of  small  drop  size  and  non-selective 
scattering  in  the  large  particle  limit,  where  the  efficiency  factor  tends  to  a  constant  value 
close  to  two. 

The  transmitted  beam  irradiance  here  is  similar  to  Rayleigh  scattering  and 

is  given  by: 


7  =  /0  exp(-<7fli?)  (2.15) 

where: 

Io=  irradiance  of  the  original  radiation 
1=  irradiance  of  the  remaining  radiation 

<Ja  -  aerosol  attenuation  coefficient  (function  of  aerosol  density) 

R=  path  length 


21 


E.  POLARIZATION 


1.  Polarization  Physics 

The  polarization  of  light  is  one  of  its  fundamental  properties,  the  others  being  its 
intensity,  wavelength  and  coherence.Christian  Huygens  was  the  first  to  suggest  that  light 
was  not  a  scalar  quantity.  This  vectorial  nature  of  light  is  called  polarization.  Light 
propagating  in  air  or  space  is  a  transverse  electromagnetic  wave.  Its  electric  field  is 
perpendicular  to  the  magnetic  field  and  they  both  are  mutually  orthogonal  to  the 
propagation  direction.  In  the  case  that  the  plane  of  vibration  of  the  electric  vector  is 
spatially  fixed,  the  wave  is  said  to  be  linearly  polarized.  Under  the  conditions  where  the 
electrical  vector  rotates  in  such  a  way  that  an  ellipse  is  created  when  viewed  from  the 
direction  of  propagation,  then  it  is  said  to  be  elliptically  polarized  (with  the  special  case 
of  circular  polarization  being  included  in  this  general  category).  The  fundamental  laws 
and  polarization  equations  describe  the  basic  parameters  of  polarization  of  the  sea  and 
sky  background  as  well  as  the  target. 

a)  Stokes  Parameters. 

Sir  G.  Stokes  found  that  the  polarization  behavior  could  be  completely 
represented  in  the  terms  of  four  measurable  quantities  known  as  the  Stokes  polarization 
parameters.  The  Stokes  polarization  parameters  for  a  plane  wave  are  [Ref.5]: 


I  =  EZS+E] 

S  p 

(2.16) 

M  —  E]  —  Ez 

(2.17) 

C  =  2EEncosS 

5  p 

(2.18) 
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(2.19) 


S  =  2EEn  sin£ 

i  p 


where: 

I  =  the  total  intensity  of  the  light 

M  =  the  intensity  of  linear  s  or  p  polarized  light 

Es=  the  normal  to  the  plane  of  propagation  components  of  the  wave 
electric  field 

Ep-  the  parallel  to  the  plane  of  propagation  components  of  the  wave 
electric  field 

C  =  the  intensity  of  linear  ±45°  polarized  light 

S  =  the  intensity  of  right  or  left  circular  polarization  contained  within  the 
beam 

8  =  the  phase  difference  between  the  s  and  p  components  of  the  E  field. 

The  first  parameter  (7 )  expresses  the  total  intensity  of  the  optical  field  and 
the  other  three  (M  ,C,S)  describe  the  polarization  state.  Stokes’  parameters  can  be 
applied  not  only  to  completely  polarized  light  but  also  to  partially  polarized  light  and 
unpolarized  light.  The  degree  of  polarization  is  the  ratio  of  the  intensities  of  the  polarized 
to  the  total  intensity: 


P  = 


_  IP0i  _(M2+C2+S2) 


2x1/2 


,0<P<1 


(2.20) 
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where: 


1  pol  =  the  intensity  of  the  sum  of  the  polarization  components 

Itol  =  the  total  intensity  of  the  beam 

The  effects  that  result  from  polarization  depend  on  the  emission  as  well  as 
the  reflection  properties  of  the  various  materials. 

b)  Emission  Polarization 

The  thermal  energy  emitted  by  various  targets  and  their  background  is 
inherently  polarized.  For  sea  targets  the  emission  characteristics  of  the  sea  and  the  sky 
background  is  of  great  importance  since  a  significant  contrast  improvement  can  be 
obtained  by  filtering  off  the  background. 

( 1 )  Sea  Emission  Polarization. 

The  mathematical  models  describing  sea  surface  polarization  are 
complex  and  are  analyzed  in  great  detail  in  Reference  6  and  10.  What  is  important  for 
this  work  is  the  fact  that  both  emitted  and  reflected  components  of  the  sea  surface 
radiance  will  show  polarization,  predominantly  in  orthogonal  directions.  For  most 
materials  there  is  only  a  small  range  of  angles  at  which  elliptical  polarization  occurs. 
Even  in  those  cases  the  phase  difference  between  orthogonal  components  is  small  enough 
(due  to  small  dielectric  constants)  that  an  approximation  of  linearly  polarized  light  at  all 
angles  is  valid  [Ref. 7]. 
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EMISS1VITY 


(2.21) 


where: 


£u  =  emissivity  of  each  facet  referred  to  the  horizontal  direction 
=  emissivity  of  each  facet  referred  to  the  vertical  direction. 

The  results  from  Figure  6  hold  true  for  all  wavelengths  less  than  or  equal 
to  10  jjm ,  excluding  the  absorption  bands  [Ref. 7].  It  should  be  noted  that  the  maximum 
values  for  polarization  are  observed  at  angles  near  the  grazing  angle. 

(2)  Sky  Emission  Polarization 

Infrared  radiation  emitted  from  the  sky  is  unpolarized  [Ref .9].  As  a 
result  there  will  be  no  improvement  in  the  contrast  measurements  when  the  sky  is  viewed 
directly.  What  should  be  mentioned  though  is  that  when  the  sky  is  viewed  towards  the 
zenith  its  emission  is  less,  due  to  the  fact  that  the  atmosphere  contains  of  water  vapor  and 
aerosols  (absorbers)  that  do  act  as  good  emitters  but  have  smaller  numbers  densities.  So 
the  sky  radiance  increases  along  with  the  zenith  angle  toward  the  horizon  to  that  of  a 
blackbody  with  atmospheric  temperature.  This  is  due  to  the  maximum  thickness  of  the 
atmosphere  encountered  at  these  angles  [Ref.  10]. 

(3)  Target  Emission  Polarization 

It  has  been  demonstrated  that  painted  surfaces  such  as  those  found 
on  ships  and  other  vehicles  display  emission  polarization.  In  Figure  7  [Ref.l  1]  a  series  of 
various  paint  measurements  is  presented  ranging  from  very  smooth  to  sand  paint 
mixtures. 
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Figure  7  -  Measured  polarized  signature  components  in  the  7.5-12  pm  band  for  selected  paint 

samples  viewed  at  45°  from  normal  [Ref.l  1]. 


The  line  for  an  ideal  surface  denotes  the  degree  of  polarization  that  would 
be  seen  from  a  perfectly  specular  reflector.  The  rougher  the  sample  the  less  is  the  degree 
of  polarization.  The  degree  of  polarization  depends  as  well  on  the  index  of  refraction  of 
the  paint  and  any  possible  degradation  due  to  the  roughness  of  the  surface. 


c)  Reflection  Polarization 
From  the  total  power  law  it  follows  that: 


a+  p+r= 1 


(2.22) 
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where: 


a-  Absorptivity 
p  =  Reflectivity 
T  =  Transmissivity 

For  a  perfectly  opaque  body,  r  =0  so: 

e  +  p  =  1  (2.23) 

The  emissivity  was  treated  in  Section  2.b.3.  The  reflectance  p  will  be 
degraded  by  the  absorption  properties  of  the  specific  body. 

( 1 )  Sea  Reflection  Polarization 

By  applying  the  Fresnel  equations  to  a  Cox-Munk  sea  surface 
slope  distribution  the  following  equations  are  derived  [Ref.8]: 


£±=1-Pl 

(2.24) 

&U  =  1  ~  Pll 

(2.25) 

e  1  +  ell 
c 

(2.26) 

2 

+  C 
H 

II 

Cl. 

(2.27) 
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e  =  \-  p 


(2.28) 


The  sea  surface  acts  as  an  opaque  body.  The  percentage  polarization  due 
to  reflection  is  given  by: 


p  =  100x  —  (2.29) 

Pl+Pu 

(2)  Sky  Reflectance 

The  sky  appears  to  be  a  poor  reflector  due  to  its  very  good 
absorptance  (because  of  the  atmospheric  water  and  aerosols). 

(3)  Target  Reflectance 

The  total  polarization  coming  from  a  thermal  target  will  depend  on 
the  balance  between  thermal  emission  and  reflection  of  incident  radiance.  The  maximum 
results  due  to  reflection  polarization  will  occur  nearer  to  grazing  incident  angles. 

2.  Polarization  Filtering  Techniques 

The  method  of  polarizing  filter  construction  has  in  the  past  been  mainly  by  using 
anisotropic  crystals  or  wire  grids.  The  wire  grid  technique  has  been  applied  in  recent 
years  by  using  lithographic  techniques.  That  is,  a  conducting  grid  is  deposited  on  a 
transparent  substrate. 

In  this  way  the  filter  will  block  the  incoming  radiation  with  the  electric  field 
parallel  to  the  orientation  of  the  grid.  That  is  because  the  wire  grid  induces  currents  that 
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reflect  the  part  of  the  incoming  wave  with  an  electric  field  parallel  to  the  wire  grid.  In 
order  for  this  process  to  take  place  the  spacing  of  the  wires  in  the  wire  grid  must  be  of  the 
order  of  the  wavelength  of  the  incident  radiation. 

For  the  infrared  region  typical  transparent  substrates  are  materials  such  as  Silicon, 
Germanium  and  KRS-5  (the  substrate  for  the  filters  actually  used  in  the  experiment) 
[Ref  .26]. 


F.  TARGET  TO  BACKGROUND  CONTRAST 

The  basic  parameter  that  is  used  by  a  system  or  an  operator  in  order  to  distinguish 
between  a  target  and  its  surroundings  is  the  contrast.  The  contrast  is  the  difference 
between  the  irradiance  received  from  a  target  pixel  and  the  irradiance  received  coming 
from  the  neighboring  background. 

There  are  several  definitions  and  types  of  contrast  used.  What  is  known  as 
inherent  contrast  [Ref.  12],  or  contrast  relative  to  the  background  is  defined  by: 


N  —  N 

_  1  T  tgt  background 


N, 


background 


(2.30) 


where: 

Ntgt  =  irradiance  coming  from  the  target 
N background  =  irradiance  coming  from  the  background. 


An  alternative  definition  of  contrast  is  presented  by  the  following  formula: 


N  -N 

£  _  tgt  backgroung 

0  ~  M  +  AT 

tgt  background 


(2.31) 


The  absolute  contrast  is  given  by: 


AN  N  lgt  N  background 


(2.32) 


When  comparing  two  images,  it  is  useful  to  quantify  the  increase  or  decrease  in 
image  contrast  improvement.  So  the  improvement  factor  is  defined  as: 


A AT, 
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(2.33) 


The  contrast  improvement  factor  for  polarization  over  the  unpolarized  case  is 
given  by: 


(2.34) 


where: 

Cpol  =  contrast  of  the  image  scene  with  polarizer 
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Cunpol  =  contrast  of  the  image  scene  without  polarizer. 

It  has  been  shown  that  the  target  to  background  contrast  outside  the  sun  glint 
corridor  can  be  improved  with  the  use  of  horizontal  polarization  filters  [Ref.  13]. 
Although  there  is  a  reduction  of  the  incoming  radiation  due  to  the  polarizing  filters  there 
was  an  improvement  in  the  contrast.  With  the  use  of  a  vertically  polarized  filter  much  of 
the  predominantly  vertical  polarized  energy  coming  from  the  background  was  rejected. 


G.  FORWARD  LOOKING  INFRARED  SYSTEMS 

A  Forward  Looking  Infrared  System  (FLIR)  is  the  infrared  analog  of  a  visible  TV 
camera  [Ref.  14].  Its  purpose  is  to  detect  and  process  the  natural  radiation  that  every 
material  body  emits.  Planck’s  laws  describe  this  kind  of  radiation  known  as  "thermal 
radiation".  The  atmospheric  transmission  windows  that  allow  the  transmission  of  thermal 
radiation  are  from  3  /zm  up  to  5  jum  (known  also  as  medium  wavelength  infrared  or 
MWIR)  and  from  8  jjm  to  12  fjm  (known  also  as  long  wavelength  infrared  or  LWIR). 

It  should  be  made  clear  that  the  FLIR  uses  different  principles  of  operation  from 
systems  that  operate  in  near  darkness  by  amplifying  the  existing  visible  (or  near  visible) 
external  radiation  due  to  moonlight,  starlight  or  any  other  such  sources.  These  systems 
are  known  as  "night  vision  devices".  Night  vision  devices  are  smaller  and  simpler  than 
the  various  FLIR  devices  but  they  do  have  certain  restrictions  such  as  the  fact  that  they  do 
not  work  satisfactorily  when  there  is  not  sufficient  light  or  in  the  presence  of  deep 
shadows.  Also  their  performance  is  inferior  to  that  of  a  FLIR  in  conditions  of  smoke, 
haze  and  dust.  The  connotation  of  a  FLIR  is  a  system  that  provides  a  visible  rendition  of 
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the  invisible  infrared  radiation  scene  [Ref.l4:p.l05].  This  technology  has  various 
applications  such  as  missile  guidance,  threat  warning,  search  and  track  and  others. 

The  purpose  of  a  FLIR  system  is  to  form  a  real  image  of  the  infrared  scene, 
identify  the  differences  in  the  imaged  radiation,  and  after  some  signal  processing 
procedure,  to  create  a  visible  representation  of  this  variation  (for  the  man  in  the  loop 
systems). 


Figure  8  -  Modular  FLIR  concept.  Common  module  systems  scan  an  array  of  IR  sensitive 
detectors  across  the  target  scene.  The  detector  outputs  are  used  to  drive  a  geometrically 
equivalent  array  of  light-emitting  diodes  to  present  a  visible  representation  to  a  television 
camera  [Ref.  14]. 
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In  early  1970’s  the  US  Department  of  defense  suggested  a  modular  approach  to 
thermal  imaging  system  design.  This  had  as  a  result  the  creation  of  the  "common 
module"  components  that  have  influenced  military  design  ever  since.  The  basic 
philosophy  of  the  common  module  system  was  to  standardize  the  basic  building  blocks  of 
the  infrared  imaging  systems.  A  first  generation  common  module  system  is 
schematically  represented  in  Figure  8  [Ref.  14]. 

These  "first  generation"  common  module  FLIRS  are  now  being  replaced  by 
systems  using  scanning  integrated  Focal  Plane  Arrays.  The  procedures  taking  place  in  an 
infrared  imaging  system  can  be  grouped  into  four  basic  categories:  detection,  image 
formation,  signal  processing  and  display. 

1.  Infrared  Detection 

There  are  two  ways  by  which  infrared  radiation  can  be  detected  by  its  interaction 
with  matter.  First  through  the  process  of  photon  exchange,  which  can  produce  heating 
and  second  by  direct  electron  emission  by  photon  absorption  (internal  photoelectric 
effect). 

In  the  photon  to  phonon  exchange  process,  the  radiation  is  absorbed  as  heat, 
resulting  in  a  rise  in  temperature.  This  increase  in  temperature  has  as  a  result  an  effect  on 
some  temperature  dependent  parameter  of  the  detector.  This  change  to  a  measurable 
physical  parameter  is  then  taken  as  a  reference  to  indicate  the  resultant  temperature 
change.  These  devices  are  known  as  square  law  detectors. 

In  the  process  of  direct  photon  absorption,  electrons  from  the  valence  band  of  a 
band  gap  semiconductor  are  elevated  to  the  conduction  band  by  direct  photon  absorption. 
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This  increase  in  the  number  of  electrons  in  the  conduction  band  is  measurable.  These 


devices  are  known  as  quantum  detectors. 

The  most  impressive  difference  between  quantum  and  square  law  detectors  is  that 
the  first  must  be  cooled  to  cryogenic  temperatures  [Ref.  14],  because  the  very  small  gaps 
required  for  infrared  photon  energies  result  in  very  high  (to  unacceptable)  numbers  of 
conduction  band  electrons  at  temperatures  close  to  ambient  and  necessitate  operation  at 
80K  or  lower.  Quantum  detectors  are  therefore  mounted  to  the  cold  surface  of  a  cooler 
inside  a  vacuum  enclosure  called  a  dewar.  However  the  sensitivity  is  greater  by  more 
than  a  factor  of  10  than  thermal  detectors. 

2.  Image  Formation 

The  electromagnetic  radiation  laws  of  reflection,  refraction  and  diffraction  for  the 
infrared  are  the  same  as  for  visible  light.  The  concepts  and  analytical  techniques  are 
directly  applicable  and  most  visible  components  have  an  infrared  counterpart 
[Ref.l4:p.l  12].  While  reflective  materials  for  the  infrared  and  visible  light  are  almost  the 
same,  refractive  elements  (such  as  visible  glasses)  are  not.  FLIR  optical  systems  also  tend 
to  be  bigger  than  visible  systems  with  the  same  field  of  view  because  of  diffraction  and 
sensitivity  differences  (depending  on  the  ratio  A/D). 

The  system  used  commonly  requires  an  optical  scanning  subsystem  to  provide 
full  focal  plane  coverage  with  limited  number  of  detector  elements  [Ref.  14].  This  is  done 
with  the  use  of  a  mechanical  scanning  system,  which  scans  the  image  across  the  detector 
in  a  prescribed  manner.  Usually  a  mirror  intersects  the  optical  beam  and  rotates  or  rocks 
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to  provide  the  required  scan.  Almost  all  scanning  systems  use  quantum  detectors.  Thus 
there  are  performance  considerations  that  require  reducing  the  background  radiation  to 
the  minimum  possible  level.  The  minimum  is  the  radiation  coming  from  the  image 
background.  In  order  to  achieve  this  the  field  of  view  of  the  detector  array  must  be 
limited  by  a  cryogenic  enclosure.  This  limiting  cryogenic  enclosure  is  called  cold 
shielding.  The  cold  shield  must  be  at  low  temperature,  or  else  it  will  add  some  part  of  the 
thermal  energy  coming  from  places  other  than  the  image  scene. 

3.  Image  Processing 

"The  function  of  the  signal  processing  of  the  signal  is  basically  to  condition  the 
detector  outputs  into  a  form  suitable  for  interface  to  the  display"  [Ref.  14].  On  the  other 
hand  the  information  processing  subsystems  function  is  to  extract  information  from  the 
signal  processing  output.  Such  functions  include  preamplification  of  the  detector  signal  in 
order  to  increase  the  level  of  the  signal  and  dc  level  subtraction  to  eliminate  unwanted 
constant  outputs  associated  with  the  background  radiation. 

In  the  common  module  systems,  each  detector  output  is  brought  out  of  the  dewar 
as  a  separate  output.  Analog  signal  processing  is  used  throughout,  one  channel  per 
detector.  Multiplexing  and  reformatting  are  done  in  the  subsystem  called  the  electro- 
optical  multiplexer.  Then  the  output  of  each  channel  is  used  to  drive  a  visible  Light 
Emitting  Diode  (LED)  in  an  array  geometrically  equivalent  to  the  detector  array.  Then 
the  diode  array  is  scanned  (a  task  performed  by  the  backside  of  the  infrared  scan  mirror 
as  is  shown  in  Figure  8)  and  is  focused  by  a  lens  (operating  in  the  visible  part  of  the 
spectrum)  to  a  TV  camera. 
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4.  Display 

The  display  is  used  as  a  converter  of  the  signal  processing  output  into  a  two 
dimensional  visible  light  presentation  of  the  selected  part  of  the  infrared  scene.  The 
display  could  be  a  standard  TV  in  the  form  of  a  cathode-ray  tube.  Care  must  be  taken  in 
this  part  as  well  as  the  previous  mentioned  components  of  the  FLIR  system,  so  that  there 
is  no  degradation  in  the  performance  of  the  system.  It  must  be  also  mentioned  that  the 
display  is  not  a  necessary  component  of  all  systems  since  there  are  applications  with  no 
man  in  the  loop.  In  those  systems  information  does  not  have  to  be  presented  in  a  form 
that  the  human  brain  is  able  to  understand  (i.e.  TV  image)  as  it  is  the  case  of  automated 
(Automatic  Target  Recognition)  systems. 
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III.  EXPERIMENT  DESCRIPTION  AND  SCENARIOS 


In  this  chapter  the  experiment  and  the  equipment  that  was  used  will  be  described. 
The  scenarios,  including  data  from  the  actual  experiment  that  in  turn  are  used  for  the 
thesis  work,  will  be  also  described.  The  EOPACE  data  were  acquired  in  April  1996  in 
San  Diego  Bay.  A  total  of  5,508  infrared  images  were  taken  of  the  R/V  POINT  SUR  in 
April  9  and  April  10,  using  an  AGA-780  Thermovision  camera  that  was  positioned  at  a 
building  at  Point  Loma.  Internal  polarizer  filters  were  used  to  get  sets  of  images  for  three 
different  polarization  cases  (horizontal,  vertical  and  unpolarized).  All  measurements 
were  taken  in  the  Long  Wavelength  Infrared  band  (8-12  fim).  Meteorological  data, 
Global  Positioning  System  (GPS)  data  and  ship  skin  temperature  data  were  recorded  by 
instruments  installed  onboard  the  ship.  In  Appendix  A  a  map  of  the  data  acquisition  area 
(San  Diego  Bay)  is  presented  as  well  as  the  GPS  ship  movement  during  the  experiment. 


A.  EXPERIMENT  EQUIPMENT 

1.  Internal  Polarizing  Filtering  with  the  AGA-780 

The  camera  that  was  used  for  the  experiment  was  the  AGA-780  Thermovision.  It 
is  a  two-channel  infrared  scanning  single  detector  system.  It  can  operate  in  both  the 
MWER  (3-8 fun)  and  LWIR  (8-14 fim)  bands  using  liquid  nitrogen  cooled  InSb  and 
HgCdTe  detectors,  respectively  [Ref.  15].  The  nominal  working  temperature  is  77  K. 
Scanning  is  accomplished  by  two  rotating  prisms  (vertical  and  horizontal  scanning  prisms 
in  Figure  9)  in  each  channel  which  scan  the  image  across  the  detector  yielding  1.1  mrad 
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geometric  resolution.  The  input  lenses  for  the  long  wave  infrared  (LWIR)  unit  used  for 
the  experiments  consist  of  anti -reflection  coated  germanium.  Its  optical  system  is  f/ 1.87, 
with  a  7x7-degree  field  of  view  [Ref.15].  Interchangeable  lenses  with  a  3.5°x3.5°  FOV 
are  also  available. 


Figure  9  -  Internal  layout  of  the  AGA-780  Thermovision  [Ref.15]. 


The  internal  filtering  is  accomplished  after  the  radiation  has  passed  through  the 
lenses.  The  polarizing  filters  (KRS-5)  are  placed  on  the  selectable  filter  wheel,  between 
the  two  components  of  the  collimation  lens,  inside  the  camera.  The  original  purpose  of 
the  filter  wheel  was  to  hold  various  accessory  spectral  filters  in  the  optical  path,  before 
the  radiation  reached  the  detector  [Ref.15].  Polarization  measurements  with  the  AGA-780 
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require  sequential  images  with  adjacent  filter  wheel  positions  inserting  vertically  and 
horizontally  polarized  filters.  The  requirement  of  turning  the  filter  wheel  in  order  to  use 
the  different  polarizer  filters  is  the  reason  for  the  sequential  nature  of  the  imaging 
procedure  during  the  experiment. 

The  two  polarizing  filters  that  were  used  each  consisted  of  an  aluminum  grid 
superimposed  on  a  9.5mm  diameter  x  3mm  thick  round  KRS-5  substrate  manufactured 
by  Graseby-Specac  (Suffolk,  England).  The  performance  specifications  provided  with  the 
filters  are  shown  in  Appendix  B.  The  different  traces  numbered  there  as  1,  2  and  3 
represent  power  transmitted  with  the  grid  vertical,  horizontal  and  with  the  two  grids  from 
the  two  filters  crossed,  relative  to  the  electronic  vector  of  the  incoming  radiation.  The 
determination  of  the  parallel  and  perpendicular  transmittances  for  the  filter  requires  the 
use  of  a  spectrometer  which  itself  includes  internal  polarizing  reflections.  So  it  must  be 
done  by  measuring  the  transmittance  two  times,  with  the  grid  in  two  orthogonal 
orientations.  These  two  measured  transmittances  are  designated  Ey  and  Eh.  The 
corresponding  transmittance  of  the  camera  is  measured  without  the  grid  in  place  as  Eo 
and  the  transmittance  for  two  crossed  grids  as  Ei.  From  these  quantities  the  grid 
transmittances  are  given  by: 
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The  manufacturer  provided  those  quantities  for  each  filter.  It  must  be  also 
mentioned  that  according  to  the  manufacturer,  for  the  traces  in  Appendix  B  the 
spectrometer  output  was  normalized  such  that  Eo=1.00  (100%)  for  all  wavelengths. 

2.  Research  Vessel  POINT  SUR 

Research  Vessel  (R/V)  POINT  SUR  is  owned  by  the  National  Science  Foundation 
(NSF).  It  is  operated  for  the  Central  California  Oceanographic  Cooperative  (CENCAL) 
by  Moss  Landing  Marine  Laboratories,  which  is  on  Monterey  Bay.  Its  equipment  consists 
of  the  necessary  navigational,  laboratory  and  mechanical  facilities  that  support  biological, 
geological,  chemical  and  physical  oceanographic  research.  For  the  purposes  of  this 
project  it  was  used  as  a  target  for  infrared  imaging  as  well  as  for  atmospheric  sampling 
functions.  In  the  Appendix  C  are  presented  the  ship’s  general  characteristics  and 
dimensions  (R/V  POINT  SUR  Cruise  Planning  Manual). 

3.  Software 

a)  CEDIP  PTRWIN 

The  analog  data  provided  by  the  scanning  detector  were  digitized  for 
display,  recording  and  measurement  by  PTRWIN  version  3.16  from  CEDIP.  This  data 
acquisition  and  analysis  software  uses  a  12-bit  A/D  converter  for  converting  the  detector 
signal  intensity  into  a  real-time  image  on  a  computer  screen  and  for  digital  image 
recording.  In  addition,  the  PTRWIN  software  was  used  as  a  visualization  tool  to  select 
images  to  be  analyzed.  Although  PTRWIN  also  offers  good  analysis  tools,  it  works 
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interactively  analyzing  only  one  image  a  time  [Ref.  16].  Following  the  process  developed 
for  a  previous  thesis  [Ref.  16],  other  software  will  be  used  for  the  analysis  and  data 
retrieved  from  the  collected  images. 

b)  Interactive  Data  Language 

IDL  (Interactive  Data  Language)  is  a  complete  computing  environment  for 
the  interactive  analysis  and  visualization  of  data.  IDL  integrates  an  array-oriented 
language  with  numerous  mathematical  analysis  and  graphical  display  techniques 
[Ref.  17].  The  image  analyses  in  this  project  were  performed  using  programs  written  in 
IDL  5.0.3  from  Research  Systems  Inc.,  Boulder,  CO.  The  programs  used  were  originally 
written  or  adapted  for  the  purpose  of  another  thesis  [Ref.  16]  and  then  modified  to  suit  the 
needs  of  this  work. 

c)  SeaRad 

SeaRad  is  a  modified  version  of  the  U.S.  Air  Force  program 
MODTRAN2,  which  predicts  the  sea  radiance  (brightness)  for  any  viewing  geometry  in 
the  spectral  range  from  52.63  cm'1  to  25000cm'1  in  addition  to  transmittance  and  path 
radiance.  The  particular  SeaRad  code  used  for  this  thesis  can  also  be  used  to  predict  the 
polarization  of  sea  radiance  [Ref.28]. 

Sea  radiance  is  modeled  by  combining  the  methods  of  geometrical  optics 
with  the  Cox-Munk  statistical  description  for  wind  driven  ocean  capillary  waves.  The 
model  is  incorporated  into  the  atmospheric  transmittance/radiance  code  MODTRAN2  to 
provide  numerical  sea  radiance  predictions.  In  the  model  each  individual  capillary  wave 
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facet  is  allowed  to  reflect  the  sky  or  sun  and  emit  thermal  radiation.  The  total  radiance 
from  the  sea  is  obtained  by  applying  the  proper  statistical  weight  to  each  facet  and 
integrating  it  over  all  facets  within  the  observer’s  field  of  view  [Ref.  18].  SeaRad  operates 
exactly  like  the  original  MODTRAN2  code  except  that  a  new  logical  parameter  (Sea 
Switch)  is  required  in  the  input  file  (tape).  Sun  glint  is  included  in  the  sea  radiance 
prediction  provided  that  the  user  has  chosen  to  execute  SeaRad  in  radiance  mode  with 
solar  scattered  radiance  included  (IEMSCT=2). 

It  should  be  mentioned  that  preliminary  comparisons  show  that  SeaRad 
agrees  to  within  a  few  degrees  Celsius  with  the  actual  sea  radiance  measurements  in  the 
mid  and  long-wave  infrared  bands.  Appendix  D  provides  an  example  of  how  SeaRad  is 
used  to  predict  ocean  radiance. 


B.  THE  EXPERIMENT 

The  experiment  was  conducted  in  San  Diego  Bay  at  Point  Loma.  The  basic 
instruments  in  the  experiment  setup  were  the  sensor  (camera  AGA-780  Thermovision) 
and  the  target  (R/V  POINT  SUR).  The  camera  was  located  at  building  15  in  NRaD,  Point 
Loma  with  coordinates  N  32°  39’  36”,  W  1 17°  14’  22".  The  ship  positioning  was  planned 
to  correspond  to  a  pattern  of  radial  bearings  with  stationary  points  spaced  at  0.5  nmi 
intervals.  At  each  stationary  point,  the  ship  should  maneuver  to  be  at  different  heading  for 
each  image  sequence  (3  polarization  cases)  taken.  Both  ship  GPS  location  and  sensor 
position  are  known,  so  the  bearing  is  computed.  From  ship  heading  and  bearing 
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information,  the  target  aspect  angle  can  be  calculated.  In  Appendix  A  the  area  of  the 
experiment  is  presented. 


C.  DATA 

1.  Images 

The  images  were  recorded  sequentially  for  the  three  polarization  cases: 
horizontal,  vertical  and  unpolarized.  An  average  of  26  images  (frames)  was  taken  for 
each  polarization  case.  Horizontal  and  vertical  cases  were  recorded  in  the  same  file  and 
each  unpolarized  case  was  recorded  in  a  different  file.  The  frames  are  composed  of  4 
interlaced  fields  with  dimensions  125x64  pixels  each.  Therefore,  each  frame  is  an  array 
of  125x256  pixels.  Image  data  taking  was  implemented  with  the  AGA-780  mounted  on  a 
remotely  controlled  pan-tilt  head  outside  the  south  looking  window  of  Building  15  on 
Point  Loma.  The  data  were  digitized  and  recorded  with  the  PTRWIN  acquisition  boards 
installed  in  the  docking  station  of  an  NEC  laptop  computer.  For  storage,  the  files  were 
transferred  to  a  second  portable  computer  and  stored  on  optical  disk.  The  original  image 
files  from  PTRWIN  are  identified  by  the  extension  *.PTW  [Ref.  16]. 

With  the  use  of  the  modified  program  BASEGEN.PRO  (written  in  IDL)  from 
Ref.  16,  140  original  image  files  were  converted  into  basic  image  files  (*.PTE).  The  basic 
files  were  named  sequentially  using  the  name  BASE  followed  by  a  number.  For  this 
thesis  data  for  the  scenarios  were  taken  from  file  BASE33.PTE  due  to  the  time  proximity 
with  radiosonde  data,  that  in  turn  were  used  for  the  construction  of  the  atmospheric 
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model  in  SeaRad.  The  structure  of  the  basic  files  for  the  three  polarization  cases  is 
presented  in  Table  4. 


Table  3  -  Original  Image  Files  (*.PTW)  Structure  [Ref.  16]. 


DATA 

OFFSET  (bytes) 

Main  Header 
(65,536  bytes) 

General  Parameters 

0 

Scanner  Parameters 

44 

Display  Parameters 

245 

Digitalization  Parameters 

375 

Comments 

563 

Setup  File  Name 

1563 

1st  Frame 
Header 
(1,536  bytes) 

Level 

0 

Time 

80 

Comments 

180 

1st  Frame 
(64,000  bytes) 

1st  Field 

0 

2nd  Field 

16,000 

3rd  Field 

32,000 

4th  Field 

48,000 

2nd  Frame 
Header 
(1,536  bytes) 

Level 

0 

Time 

80 

Comments 

180 

2nd  Frame 
(64,000  bytes) 

1st  Field 

0 

2nd  Field 

16,000 

3rd  Field 

32,000 

4th  Field 

48,000 

The  Frame  Structure  above  is  repeated  until  the  last  frame  in  the  image  sequence. 


46 


Table  4  -  Basic  Files  (*.PTE)  structure.  [Ref.16]. 


Variable 

Type 

Name  of  original  (*.PTW)  file 

string  (20) 

Information  vector 

integer  (1,13)' 

Frame  time  :  horizontal  polarization 

integer  (4,  #  frames  hor.  case) 

Frame  time  :  vertical  polarization 

integer  (4,  #  frames  ver.  case) 

Frame  time :  unpolarized 

integer  (4,  #  frames  unp.  case) 

Image  set  (frames) :  horizontal  polarization 

integer  (125, 256,  #  frames  hor.  case) 

Image  set  (frames) :  vertical  polarization 

integer  (125, 256,  #  frames  ver.  case) 

Image  set  (frames) :  unpolarized 

integer  (125,  256,  #  frames  unp.  case) 

Results :  horizontal  polarization 

float  (100,  #  frames  hor.  case) 

Results :  vertical  polarization 

float  (100,  #  frames  ver.  case) 

Results :  unpolarized 

float  (100,  #  frames  unp.  case) 

2.  Meteorological  Data 

The  NPS  Boundary  Layer  Meteorology  Group  collected  meteorological  data  on 
two  meteorological  towers  (METOC  1  and  METOC  2)  installed  onboard  the  ship.  The 
data  were  recorded  every  30-seconds.  Ship  GPS  data  are  included  in  the  METOC  files. 
Also  radiosonde  data  were  collected  every  four  to  five  hours.  Appendix  E  describes  and 
presents  the  format  used  to  record  the  data  collected  on  the  ship  as  well  as  the  radiosonde 
data.  The  two  meteorological  towers  were  located  onboard  the  ship  in  the  following 
configuration:  METOC  1  tower  was  above  the  ship  bridge  (approximately  38  ft  above  sea 


1  In  DDL,  the  first  matrix  index  is  column  number  and  the  second  is  line  number. 
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surface)  and  METOC  2  tower  was  at  the  ship  bow  (approximately  24  ft  above  sea 
surface). 


3.  Ship  Skin  Temperature  Data 

The  skin  temperature  of  the  ship  was  measured  with  16  thermistors  installed  by 
NPS-NACIT  at  different  points  on  the  ship.  Measurements  were  taken  and  recorded 
every  20  seconds.  In  Appendix  C  is  presented  the  position  for  each  thermistor  onboard 
the  ship  and  the  format  used  to  record  the  data. 


D.  SCENARIOS 

For  the  simulation  two  scenarios  were  designed.  The  target  used  was  an 
equivalent  parallelepiped  with  the  dimensions  (length,  width  and  height)  of  R/V  POINT 
SUR.  The  sensor  simulated  has  the  generic  characteristics  of  a  Common  Module  Forward 
Looking  Infrared  sensor  (FLER).  The  characteristics  are  given  in  chapter  IV. 

For  the  scenarios  it  was  assumed  that  the  FLIR  was  mounted  on  an  aircraft.  The 
aircraft  approached  the  target  in  level  flight  at  two  different  heights:  100m  and  1000m.  A 
schematic  representation  of  the  scenarios  is  presented  in  Figure  10.  For  the  input 
characteristics  to  the  atmospheric  model  used  in  SeaRad  the  data  from  the  meteorological 
measurements  (basically  the  radiosonde)  were  used.  The  temperature  of  the  hot  spot  of 
target  (R/V  POINT  SUR)  was  obtained  from  image  file  BASE33.PTE  with  the  assistance 
of  some  programs  in  IDL.  Those  programs  were  initially  written  for  the  purposes  of 
Ref.  16  and  then  modified  to  suit  the  needs  of  this  thesis. 
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Figure  10  -  Schematic  representation  of  the  scenarios. 
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IV.  SYSTEM  PERFORMANCE  FIGURES  OF  MERIT  AND  DETECTION 

CRITERIA 

In  this  chapter  the  FLIR  system  performance  figures  of  merit,  the  apparent 
temperature  difference  and  the  various  detection  criteria  will  be  presented. 


A.  FLIR  PERFORMANCE  MEASURES. 

1.  Noise  Equivalent  Temperature  Difference 

Noise  Equivalent  Temperature  Difference  (NET  or  NETD)  is  a  useful  expression 
for  the  sensitivity  of  a  thermal  system.  It  is  the  temperature  difference,  referenced  to  300 
K,  between  a  large  target  (a  square  that  is  at  least  ten  times  the  detector  angular  subtense) 
and  the  background,  which  is  required  to  produce  a  peak  signal-to-rms-  noise  of  one  at  a 
particular  point  in  the  signal  processing  chain  [Ref. 2]. 

The  NET  is  an  excellent  diagnostic  tool  for  production  testing  to  verify 
sensitivity.  It  is,  none  the  less,  a  poor  system  to  system  comparison  parameter  and  should 
be  used  cautiously  when  comparing  systems  built  to  different  designs.  That  is  because 
NET  is  a  function  of  spectral  responsivity  and  the  noise  equivalent  bandwidth  [Ref.  12]. 

The  form  of  NET  is  given  in  Ref.  19  as: 


(A/„)1/24(/  i#y 


NET  =  - 


JlT0TAAdV2K'2 


[  ^-D*xdZ 

JAA  TIT  *■ 


Iaa  37- 


(4.1) 
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where: 


Lx  =  spectral  radiance  (W  cm'2  sr'1  f£m~x ) 

/  /#=  system  effective  /  /# 

T0  =  optics  transmission 

ta=  atmospheric  transmission  over  the  laboratory  path 
used  to  make  measurement 
Ad  =  detector  area  (cm2) 

Ns  =  number  of  detectors  in  series 
A/l  =  spectral  band  (fJm) 

dL 

— —  =  derivative  with  respect  to  temperature  of  the 
dT 

spectral  radiance 

D\  -  spectral  detectivity  (cm  Hz1/2  W'1) 

A/„=  electronic  noise  bandpass 

A  drawback  in  the  use  of  NET  as  a  figure  of  merit  is  the  fact  that  it  does  not 
include  display  or  observer  response  characteristics.  So  NET  does  not  relate  directly  to 
human  performance. 
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a)  Minimum  Detectable  Temperature  Difference 
Minimum  Detectable  Temperature  Difference  (MDT  or  MDTD)  is  a 
laboratory  measurement  of  a  thermal  system’s  sensitivity,  which  includes  a  human 
operator  (unlike  NET). 

MDT  is  the  apparent  temperature  difference  between  a  square  (or  circular) 
target  and  a  uniform  background,  required  for  the  detection  of  the  square  (or  circular) 
target  by  a  trained  observer  viewing  through  a  FLIR.  There  are  no  time  constraints  in  the 
observation  process  and  the  location  of  the  target  is  approximately  known  to  the  observer 
[Ref. 2].  MDT  is  measured  for  a  set  of  targets  of  increasing  area  and  is  usually  plotted  as 
AT  versus  l/(target  area)172  as  in  Figure  11.  It  can  be  observed  from  Figure  11  that  the 
curve  is  not  asymptotic  for  small  targets  since  an  increase  in  temperature  can  make  even 
a  very  small  target  detectable,  even  when  the  angular  subtense  is  smaller  than  the 
instantaneous  field  of  view  (IFOV)  of  the  detector. 

Shumaker,  Wood  and  Thacker  give  a  formula  for  MDT  as  [Ref.21]: 


MDT(y)  = 


0 SNRT)(NET)(Qt  +  r 2  )(AxAy)172 


Q7 


K 


-il/2 


^(r?+rZ+nT)teFrNosNss 
4 


(4.2) 


where: 

SNRT  =  signal  to  noise  threshold 

NET  =  noise  equivalent  temperature  difference 

Q.T  =  target  angular  subtense  (mrad2) 
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rg  =  system  resolution  (mrad) 

Ax  =  in  scan  detector  subtense  (mrad) 

Ay  =  cross  scan  detector  subtense  (mrad) 

rB  =  back  end  resolution  (mrad) 

te  =  eye  integration  time  in  seconds 

Fr  =  frame  rate  in  Hz 

Nos  =  overscan  ratio 

Nss  =  serial  scan  ratio 


Figure  11  -  Typical  plot  of  MDT  showing  relative  sizes  of  the  MDT  pattern  [Ref.2:p.81]. 
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b)  Minimum  Resolvable  Temperature  Difference 
Following  previous  work  on  television  systems  Lloyd  and  Sendall 
[Ref.20]  developed  the  thermal  imaging  system  performance  measure,  Minimum 
Resolvable  Temperature  Difference  (MRT  or  MRTD). 

MRT  is  a  noise  limited  threshold  measure  of  system  spatial  resolution  and 
thermal  sensitivity,  in  which  the  characteristics  of  the  human  observer  are  included  (thus 
making  it  a  subjective  measurement).  It  is  defined  as  the  temperature  difference  between 
a  uniform  background  and  the  bars  of  a  four  bar  pattern,  each  bar  having  a  7:1  aspect 
ratio  (so  the  overall  pattern  will  be  a  square),  which  is  required  by  a  trained  observer  to 
just  resolve  all  four  bars  when  viewing  the  pattern  through  the  imager  [Ref.2].  A  example 
of  an  MRT  plot  is  displayed  in  Figure  12. 

As  can  be  seen  from  Figure  12,  the  temperature  difference  that  is  needed 
to  resolve  the  four  bars  increases  as  the  bars  get  smaller.  The  MRT  curve  is  finally 
asymptotic  at  a  spatial  frequency  near  Ax-1  (unlike  MDT),  since  even  very  large 
temperature  differences  will  not  allow  resolution  of  the  four  bar  patterns  where  one  cycle 
of  the  bar  pattern  is  smaller  than  a  detector  angular  subtense. 

Shumaker  provides  a  formula  for  the  MRT  [Ref.21  :p.8-52]: 


„nrr,  s  2  xSNRTxNETxp1'2 
MRT  (v)  = - -x 


MTFAv) 


v2AxAy 


nl/2 


1-1/2 


(4.3) 


where: 


SNRT  =  signal  to  noise  threshold 

NET  =  noise  equivalent  temperature  difference 

px  =  noise  filter  factor 

v  -  spatial  frequency  for  which  MRT  is  desired 
Ax  =  in  scan  detector  subtense  (mrad) 

Ay  =  cross  scan  detector  subtense  (mrad) 

MTF  =  modulation  transfer  function 
L=  length  to  width  ratio  of  the  bar  (7) 
te  -  eye  integration  time  in  seconds 
Fr  =  frame  rate  in  Hz 
Nos  =  overscan  ratio 
Nss  =  serial  scan  ratio 
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Figure  12  -  A  typical  MRT  plot  indicating  the  sizes  of  the  MRT  four  bar  pattern  at  various 

spatial  frequencies  [Ref.2:p.82]. 


In  Ref. 21  is  also  provided  another  formula  for  the  Noise  Equivalent 
Temperature  Difference,  which  can  be  derived  from  the  formula  that  was  presented 
earlier: 


10(FOVxFOVFNosNss) 


(7tNDnJ'2DAxAyD"nsc^T0 


where: 

FOVx  =  in  scan  field  of  view  in  mrad 
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FOVy  =  cross  scan  field  of  view  in  mrad 

Fr  =  frame  rate  in  Hz 

Ngs  =  overscan  ratio 

Nss  =  serial  scan  ratio 

Nd=  number  of  detectors 

nsc  =  scan  efficiency 

Ax  =  in  scan  detector  subtense  (mrad) 

Ay  =  cross  scan  detector  subtense  (mrad) 

D  =  aperture  diameter  in  meters 

D**  =  band  average  detectivity  in  cm  Hz1/2/W 

dN  2 

- =  derivative  of  Planck’s  equation  in  W/cm  K  Sr 

ar 

T0=  transmission  of  the  optics 

B.  APPARENT  TEMPERATURE  DIFFERENCE  COMPUTATION 

The  temperature  difference  AT  between  a  target  and  its  background  is  defined  for 

zero  range  A T0  (no  attenuation).  When  the  energy  from  the  target  and  its  background  is 
modified  by  the  atmosphere  in  the  path  between  the  target  and  the  sensor,  then  it  is 
known  as  the  apparent  temperature  difference  ( A Tapp ). 
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There  are  two  basic  ways  to  compute  the  apparent  temperature  difference.  One  is 
to  assume  that  the  contrast  temperature  difference  is  attenuated  in  the  same  fashion  as  the 
radiance: 

Arw=A7>-'s  (4.5) 


where: 

A Tapp  =  apparent  temperature  difference  at  range  R 
A T0=  zero  range  temperature  difference 
H  =  extinction  coefficient 
R  =  range 

The  use  of  equation  4.5  for  the  computation  of  the  apparent  temperature 
difference  is  a  very  simple  method,  frequently  used,  and  its  results  are  not  very  accurate. 

Another  method  is  to  convert  the  temperature  difference  into  in-band  background 
radiance  and  target  radiance,  compute  the  corrections  for  attenuation  and  path  radiance 
by  using  a  specific  code  (i.e.  LOWTRAN,  MODTRAN,  SEARAD)  and  reconvert  to  an 
apparent  A Tapp  as  a  function  of  range.  This  method  is  much  more  accurate  than  the 

previous  one  and  will  be  used  for  this  thesis  along  with  the  results  from  the  SEARAD 
propagation  code. 

For  the  conversion  of  the  temperature  difference  into  in-band  background 
radiance  and  target  radiance  a  program  written  in  MATLAB  was  used  (rad.m).  The  code 
is  provided  in  Appendix  F.  This  program  performs  a  numerical  integration  of  Planck’s 
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law  over  the  bandwidth  and  provides  lookup  tables  that  associate  temperature  with 
radiance. 

The  target’s  temperature  was  obtained  from  the  file  BASE33.PTE  with  the 
assistance  of  the  modified  program  Hotspot.pro  (written  in  IDL).  The  modified 
Hotspot.pro  located  the  hot  spot  of  the  target  and  provided  the  temperature  in  degrees 
Kelvin.  Then  this  temperature  was  converted  to  radiance  with  the  program  rad.m. 

A  model  of  the  specific  atmospheric  conditions  from  the  radiosonde  data  from  the 
experiment  (Appendix  E)  was  inserted  in  SEARAD  (Appendix  D)  as  well  as  the 
geometry  for  the  two  scenarios.  The  above  data  provided  a  set  of  polarized  and 
unpolarized  values  for  the  transmittance,  path  radiance  and  background  radiance.  The 
transmittance  multiplied  by  the  ship  radiance  and  with  the  addition  of  the  path  radiance 
provided  the  total  ship  and  path  radiance  at  the  sensor.  This  radiance  was  then  converted 
to  temperature  with  the  use  of  the  Planck  formula  for  the  in-band  radiance.  The 
calculation  was  performed  by  the  rad.m  program  (Appendix  F). 

The  full  range  background  radiance  for  the  specific  geometry  (thus  a  specific 
range)  was  also  converted  to  temperature.  This  temperature  was  then  subtracted  from  the 
total  ship  and  path  radiance,  providing  the  apparent  temperature  difference  A Tapp  for  the 

specific  range. 

For  the  same  calculations  with  polarization  the  values  of  the  full  range 
background  radiance  and  the  ship  and  path  radiance  were  multiplied  by  a  factor  of  0.7 
and  0.5  respectively  to  compensate  for  the  attenuation  of  the  energy  passing  through  the 
polarizing  filter.  The  analytical  calculations  of  the  above  factors  of  0.7  and  0.5  are 
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presented  in  Appendix  B.  These  numbers  represent  the  amount  of  in-band  energy  that 
passes  through  the  polarizer  filters  for  the  cases  of  horizontally  polarized  and  unpolarized 
waves. 


C.  DETECTION  CRITERIA  AND  PROJECTED  AREA  OF  THE  SHIP 

The  initial  task  in  every  system  is  to  search  the  field  of  view  to  locate  the  target. 
The  search  pattern  varies  from  system  to  system  and  also  with  the  observer  training  and 
background.  The  following  step  is  finding  the  target.  DAgostino  and  Moulton  introduced 
the  "minimum  findable  temperature  difference"  (MFTD)  [Ref. 23].  "In  this  process  the 
target  might  be  anywhere  in  the  field  of  view.  When  the  target  is  located,  then  attention  is 
focused  on  a  particular  area  of  the  scene.  This  process  is  more  complex  than  the 
minimum  detectable  temperature  (MDT).  The  MDT  is  the  ability  to  detect  a  blob  in  noise 
with  the  target  in  the  center  of  the  field  of  view.  Both  the  MDT  and  MFTD  are  functions 
of  the  target  angular  subtense"  [Ref.  12]. 

After  finding  the  target  a  specific  set  of  criteria  have  to  be  applied  to  provide  the 
required  amount  of  detail  for  detection,  recognition  and  identification  (detection, 
recognition  and  identification  are  just  three  discrimination  levels  in  a  continuum).  The 
definition  of  each  level  depends  on  the  author  (since  each  author  has  tried  to  describe  his 
impression  of  discrimination  )  [Ref.  12].  Typically  used  levels  of  target  discrimination  are 
listed  in  Figure  13. 

Because  of  the  difficulty  of  the  description  of  a  target  mathematically,  for 
measurement  of  the  capability  of  the  FLIR  aided  eye  to  detect  (in  a  broad  sense)  a  target, 
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two  types  of  simple  targets  are  used:  isolated  rectangles  and  periodic  bar  patterns.  The 
first  step  is  the  detection  of  isolated  rectangles  and  the  second  is  the  resolution  of  bars  in 
a  bar  pattern. 


TASK 

DESCRIPTION 

EXAMLES 

Pure  Detection 

Something-Nothing  Decision 

Object  vs.  Noise 

Discrimination 

Detection 

Potential  Target  vs.  Background 
Decision 

Target  vs.  Background 

Orientation 

Determining  the  Direction  of  the 
Long  Axis  of  the  Target 

Side  View  vs.  Front  View 

Classification 

Distinguishing  Potential  Targets 
from  Similarly  Sized  Non-Targets 

Warship  vs.  Merchant  Ship 
Tracked  vs.  Wheeled  Vehicle 

Recognition 

Distinguishing  the  Target  from 
Similarly  Classed  Military  Objects 

Destroyer  vs.  Cruiser 

Tank  vs.  APC 

Identification 

Determining  the  Specific  Military 
Designation 

Kashin  vs.  DD963  Destroyer 

Unit 

Identification 

Determining  the  Identity  of  the 
Specific  Object 

USS  Kennedy  vs.  USS  Forrestal 

Figure  13  -  Discrimination  levels  [Ref.21:p.2-3]. 


The  first  task  is  described  as  "aperiodic  detection"  and  is  characterized  by  the 
FLIR  parameter  Minimum  Detectable  Temperature  (MDT).  The  second  task  is  described 
as  periodic  detection  and  is  characterized  by  the  FLIR  parameter  Minimum  Resolvable 
Temperature  (MRT). 

This  simplified  analysis  has  proven  widely  usable  as  an  analytical  approach  and  it 
is  customary  to  describe  all  visual  tasks  by  equating  them  to  one  of  the  two  problems 
above  for  which  mathematical  solutions  have  been  developed.  Generally  pure  detection  is 
equated  to  detection  of  an  isolated  rectangle,  with  area  equal  to  the  target  area,  and  whose 
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temperature  difference  with  the  background  is  set  equal  to  the  average  temperature 
difference  of  the  real  target  and  its  background. 

For  the  other  discrimination  tasks,  the  targets  are  related  with  bar  chart  patterns 
whose  temperature  difference  between  the  hot  and  cold  bars  is  equal  to  the  average 
temperature  difference  of  the  real  target  and  its  background  and  whose  spatial  frequency 
is  determined  by  the  target’s  critical  dimension  and  the  difficulty  of  the  visual  task 
attempted  [Ref.21].  This  procedure  is  illustrated  in  Figure  14.  In  the  case  where  a  specific 
feature  must  be  detected  for  the  discrimination  task,  then  its  area  alone  may  be  equated  to 
an  idealized  target. 


Ffwnupnrv  -  1  _  N>  Where  Ni  is  number  of  resolvable  cycles 
y  a”  dc  required  to  perform  a  discrimination 

(Nj  increases  with  increasing  task  difficulty) 


Figure  14  -  Target  task  transform  [Ref.21  :p.2-5]. 
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1.  Visual  Discrimination  Criteria 

There  have  been  only  a  few  experiments  designed  to  relate  the  resolution  of  bar 
charts  to  levels  of  visual  discrimination.  The  classic  experiment  was  conducted  in  1957 
by  Johnson  [Ref.24]  for  application  to  image  intensifier  systems.  The  results  from  that 
experiment  are  presented  in  Table  5  and  are  known  as  the  "Johnson  Criteria"  [Ref .21]. 


Table  5  -  Johnson  criteria 


TAR(jET  (OkOAbSlbE}  RESOLUTION  PER  MIN.  DIMENSION  (LINE  PAIRS) 


DET. 

ORIENT. 

RECOG. 

IDENT. 

■ilKIlli 

1.25 

4.5 

8.0 

M-48Tank 

0.7 

1.2 

3.5 

7.0 

Stalin  Tank 

0.75 

1.2 

3.3 

6.0 

Centurion  Tank 

0.75 

1.2 

3.5 

6.0 

Half-Track 

1.0 

1.50 

4.0 

5.0 

Jeep 

1.2 

1.50 

4.5 

5.5 

Command  Car 

1.2 

1.5 

4.3 

5.5 

Soldier  (Standing) 

1.5 

1.8 

3.8 

8.0 

105  Howitzer 

1.0 

1.5 

4.8 

6.0 

AVERAGE 

1.0  ±0.25 

1.4  ±  0.35 

4.0  ±  0.8 

6.4  ±  1.5 

Later  Moser  concluded  that  the  Johnson  Criteria  could  not  be  simply  extrapolated 
to  ship  targets  [Ref. 25].  Following  a  study  with  ship  silhouettes  he  concluded  that  the 
minimum  ship  dimension  was  of  little  importance,  but  the  classification  and  identification 
of  ships  could  be  related  to  the  number  of  pixels  within  the  area  of  the  ship  image  and 
determined  the  required  numbers  as  1,  66  and  400  pixels  at  the  image.  These  criteria 
correspond  to  the  number  of  resels  (resolution  elements)  per  average  ship  dimension  as  a 
measure  of  visual  discrimination. 
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O'Neill  used  a  television  system  and  experimented  using  a  ship  target  158  meters 
long  and  14  meters  high  with  an  experienced  observer.  Subsequent  research  at  the  Naval 
Air  Development  Center  showed  that  the  number  of  resolvable  elements  characterizing  a 
task  is  dependent  upon  the  signal  to  noise  ratio  in  the  imagery.  The  less  the  noise  the 
fewer  the  resolvable  elements  required. 

There  do  not  exist  any  criteria  for  visual  discrimination  that  hold  true  for  a  wide 
spectrum  of  targets  or  conditions.  So  the  best  way  for  the  determination  of  the 
appropriate  criteria  for  a  specific  target,  is  to  perform  a  simple  psychophysical 
experiment.  When  no  other  information  is  available  on  which  to  establish  a  bar  chart 
equivalent  the  data  in  Table  6  can  be  used  as  a  first  approach  to  determining  the 
resolution  criterion  for  the  desired  task 


2.  Critical  Dimension  and  Target  Aspect 

For  the  establishment  of  the  bar  chart  equivalent  of  a  specific  task,  the  target 
critical  dimension  must  be  established.  For  this  work  the  Moser  "average  target 
dimension",  which  is  the  square  root  of  the  target-projected  area,  will  be  used  as  the 
target  critical  dimension.  This  way  of  calculating  the  critical  dimension  was  found  by 
Moser  to  lead  to  consistent  predictions  even  for  very  elongated  targets. 
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Table  6  -  Summary  of  target  transform  data  [Ref.21]. 


ussnmm 


RESOLVED  LINE  PAIRS  (CYCLES)  PER  CRITICAL  DIMENSION 


Truck 

PET. 

0.90 

ORIENT. 

1.25 

RECOG. 

4.5 

IDENT. 

8.0 

M-48Tank 

0.75 

1.2 

3.5 

7.0 

Stalin  Tank 

0.75 

1.2 

3.3 

6.0 

Centurion  Tank 

0.75 

1.2 

3.5 

6.0 

Half-Track 

1.0 

1.50 

4.0 

5.0 

Jeep 

1.2 

1.50 

4.5 

5.5 

Command  Car 

1.2 

1.5 

4.3 

5.5 

Soldier  (Standing) 

1.5 

1.8 

3.8 

8.0 

105  Howitzer 

1.0 

1.5 

4.8 

6.0 

AVERAGE 

1.0  ±0.25 

1.4  ±  0.35 

4.0±0.8 

6.411.5 

Sfflt>  TARGETS' 

O'NEILL  DATA 


RESOLVED  LINE  PAIRS  (CYCLES) 

SION 


mm  ii  in i 

Detect  Object  on  Horizon  Sky 

3 

Recognize  a  Vessel 

5 

Recognize  Ship  Structure 

11 

Recognize  Ship  Type 

16 

Classify  King  Posts 

16 

Discern  Radar  Detail 

33 

Detea  40mm  Gun  Barrel 

55 

WM33SW; 


DISCRIMINATION  TASK 
Detect  a  Ship 

Classify  Ship  as  Combatant 
Recognize  Ship  Type _ 


RESOLVED  LINE  PAIRS  (CYCLES) 
PER  CRITICAL  DIMENSION 

(Aperiodic  treatment) 

4 

10 


Since  the  critical  dimension  will  be  calculated  as  the  square  root  of  the  exposed 
target  area,  the  aspect  of  the  target  will  play  a  major  role.  The  target  aspect  affects  the 
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apparent  size  and  shape  of  the  target.  The  area  of  the  target  image  as  a  function  of  its 
dimensions  and  the  azimuth  and  elevation  aspect  angles  is  given  by: 

AT  =  lxhxcos6xcos(p  +  wxhxcos0xsin<p  +  lxwxsin0  (4.6) 

where: 

At=  projected  target  area 

/  =  actual  target  length  in  meters 

w=  actual  target  width  in  meters 

h  =  actual  unobstructed  target  height  in  meters 

6=  elevation  angle 

<p-  azimuth  angle 

The  critical  dimension  of  the  target  then  will  be: 


Dc  — 


(4.7) 


The  orientation  of  the  target  for  the  calculation  of  the  projected  target  area  (with 
the  dimensions  of  R/V  POINT  SUR)  due  to  the  viewing  angle  is  displayed  in  Figure  15 
along  the  lines  of  Shumaker  [Ref.  21]. 
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D.  MRT  COMPUTATION  AND  CALCULATION  VS  RANGE  WITH 

VARIOUS  CRITERIA 

Minimum  Resolvable  Temperature  Difference  (MRT)  is  the  most  used  and  most 
useful  FLIR  specification  parameter.  It  gives  the  Temperature  Difference  AT  between 
hot  and  cold  bars  of  a  standard  4  bar  (7:1  aspect  ratio)  chart  required  to  make  the  bars  just 
resolvable,  as  a  function  of  the  spatial  frequency  of  the  bars  [Ref.21].  The  object  space 
spatial  frequency  fx  is  the  inverse  of  the  target  angular  subtense  and  is  usually  measured 
in  cycles  per  milliradian: 


/> 


1000 


2d 


(4.8) 


where: 

fx=  spatial  frequency  in  cycles  per  milliradian 
R  =  distance  from  the  infrared  system  entrance  aperture  to  the  target 
2d  =  the  spatial  extent  of  one  cycle  =  two  times  the  bar  width 

The  schematic  representation  of  the  above  is  displayed  in  Figure  16. 
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Figure  16  -  Definition  of  spatial  frequency.  MRT  targets  have  a  7:1  aspect  ratio  [Ref. 2], 

In  order  to  plot  the  target  spatial  frequency  as  a  function  of  range,  which  will  be 
needed  in  the  plotting  of  the  MRT  vs  range,  the  fundamental  spatial  frequency 
characterizing  the  ship  (basically  twice  its  critical  dimension)  is  expressed  as  a  function 
of  range  according  to  Figure  17.  This  procedure  is  in  accordance  with  that  described  in 
Shumaker  [Ref.21].  In  the  case  that  the  task  is  pure  detection  then  N=l.  Otherwise  a 
value  from  the  bar  chart  equivalent  has  to  be  taken. 

For  the  calculation  of  the  MRT  a  FLIR  system  with  the  characteristics  displayed 
in  Table  7  was  assumed.  The  values  displayed  in  Table  7  represent  typical  values  for  a 
generic  Common  Module  FLIR. 
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N  =  Number  of  Bars  in  the  Bar  Chart 
R  -  Slant  Range 
TAS  -  Target  Angular  Subtense 
&c  -  Critical  Dimension 
X  -  Physical  Dimension  to  be  Resolved 
*c  =  Spatial  Frequency  of  Bar 


Figure  17  -  Target  spatial  frequency  vs  range. 

Some  of  these  values  were  obtained  from  Shumaker  [Ref.21],  such  as  the  Field  of 
View  values,  the  serial  scan  ratio,  the  over-scan  ratio,  the  dwell  time  and  the  number  of 
multiplexer  samples  per  detector  dwell  time.  Others  were  taken  from  the  original  data 
from  the  experiment,  such  as  the  target’s  dimensions.  Some  were  chosen  by  assuming  that 
these  values  would  be  close  to  a  real  sensor,  since  these  values  are  not  made  publicly 
available  by  the  companies  that  build  FLIR  systems.  These  values  are  the  signal  to  noise 
ratio,  the  in  and  cross  scan  detector  subtense,  the  aperture  diameter,  the  focal  length,  the 
display  spot  size  and  the  amplitude  of  residual  image  motion.  Finally  some  of  the 
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tabulated  values  are  standard  or  were  chosen  in  accordance  to  the  specific  detection 
criteria.  These  values  are  the  number  of  bars  in  the  bar  chart,  the  specific  detectivity,  the 
average  wavelength,  the  eye  integration  time  and  the  length-to-width  ratio  of  the  bars  in 
the  bar  chart. 

From  these  parameters  the  MRT  and  NET  were  computed  using  Equations  4.3 
and  4.4,  and  expressed  as  functions  of  target  range,  using  Equation  4.8  and  the  procedure 
of  Figure  17,  for  the  values  of  N  appropriate  to  detection  (N=2)  recognition  (N=8)  and 
identification  (N=13)  using  the  Johnson  Criteria.  The  dimensions  of  R/V  POINT  SUR 
were  used  for  the  calculation  of  the  target  characteristics. 

With  the  assistance  of  a  series  of  IDL  programs  the  ship’s  hot  spot  temperature 
along  with  the  sea  surface  temperature  was  obtained  from  the  data  of  the  experiment. 

With  the  known  hot  spot  ship  temperature  (300.26  K,  the  average  of  thermistors 
11  and  16)  and  sea  temperature  (289.95  K,  sea  surface  temperature),  the  zero  range 
Temperature  Difference  (A T0  =  Ttgt  - Thckgd  )  was  converted  using  Planck’s  Law  (with  the 
rad.m  program)  into  a  Target-Background  Radiance  Difference  with  horizontal  polarizer 
(A Lh)  and  with  no  polarizer  (A L0).  These  radiance  differences  were  propagated  to 

arbitrary  range  using  SEARAD  to  provide  the  polarized  sea  surface  radiance  and  the  path 
attenuation  and  path  radiance  for  both  the  horizontally  polarized  and  unpolarized  cases. 
Reconverted  to  temperature  difference  these  are  tabulated  in  an  EXCEL  spread  sheet  (an 
example  of  one  of  spreadsheets  used  is  in  Appendix  J).  A  schematic  representation  of  the 
procedure  is  displayed  in  Figure  18. 
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Lship(R)  .  ►  Tapp.  ship(R) 


Tsea(O) - Lsea(O) 

At  zero  range 


Figure  1 8  -  Schematic  representation  of  the  procedure  followed  to  derive  the  Apparent 

Temperature  Difference  ATapp 

With  these  characteristics  of  the  system  and  the  target  and  by  taking  into 
consideration  the  target  aspect,  a  series  of  MRT  curves  were  plotted  with  the  assistance 
of  the  mrt.m  program  listed  in  Appendix  H.  The  following  plots  were  calculated  for  a 
height  of  the  observer  (sensor)  of  1000m.  N  is  the  number  of  bars  in  the  bar  chart  that  in 
this  case  is  also  the  number  of  resolution  elements  (resels,  defined  as  "the  smallest  region 
in  the  object  space  whose  dimension  is  equal  to  the  resolution  in  that  dimension" 
[Ref.  12])  required  on  the  target  dimensions. 
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Table  7  -  FLIR  system  and  target  characteristics. 


Value 

Definition 

6.86  degrees 

Field  of  View  in  the  x  direction 

FOVv 

5.16  degrees 

Field  of  View  in  the  y  direction 

2.5 

Signal  to  noise  ratio 

0.1 

In-scan  detector  subtense  in  mrad 

dy 

0.1 

Cross-scan  detector  subtense  in  mrad 

L 

7 

Length -to- width  ratio  of  the  bar  7:1 

te 

0.1 

Eye  integration  time  in  seconds 

Nss 

1 

Serial  scan  ratio 

Nos 

1 

Over-scan  ratio 

Fr 

30 

Frame  rate  in  Hz 

D 

0.1 

Aperture  diameter  in  meters 

A 

Approximation  of  the  radiation  wavelength 

1 

41.5 

Length  of  the  ship  in  meters 

w 

9.75 

Width  of  the  ship  in  meters 

h 

8.8 

Height  of  the  ship  in  meters 

td 

6*10^ 

Dwell  time  in  seconds 

7 

Frequency  at  which  the  response  is  down  3db 
from  its  low  frequency  response 

Dfp 

8.5*  106 

Half  power  frequency  response  of  the 
preamplifier  in  HZ 

Ns 

2 

Number  of  multiplexer  samples  per  detector 
dwell  time 

HD 

0.2 

Height  of  the  display  in  meters 

dss 

Display  spot  size  in  meters 

Dv 

0.5 

Display  viewing  distance  in  meters 

AR 

6.015 

Amplitude  of  residual  image  motion  in  mrad 

to 

0.7 

Transmission  of  the  optics 

Nd 

180 

Number  of  detectors 

Nsc 

0.75 

Scan  efficiency 

Ds 

Specific  detectivity  in 
[cm  Hz1/2  W'1] 

f 

0.2 

Focal  length  in  meters 

N 

Depends  on  the 
criteria 

Number  of  bars  in  the  bar  chart 
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0  2  4  6  8  10  12  14  16  18 

Range  in  meters  x  1Q4 


N=  required  number  of  resels  on  critical  dimension  for  detection 
Figure  19  -  MRT  vs  Range  for  Detection  N=2  (Johnson). 
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25 


MRT  vs  Range  for  Recognition  N=8  (Johnson) 


N=  required  number  of  resels  on  critical  dimension  for  recognition 
Figure  20  -  MRT  vs  Range  for  Recognition  N=8  (Johnson). 
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emperature  difference  in°C 


MRT  vs  Range  for  Identification  N=13  (Johnson) 


0  0.5  1  1.5  2  2.5 

Range  in  meters  x  104 

N=  required  number  of  resels  on  critical  dimension  for  identification 
Figure  21  -  MRT  vs  Range  for  Identification  N=13  (Johnson). 
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Range  in  meters 


N=  required  number  of  resels  on  critical  dimension  for  recognition 
Figure  22  -  MRT  vs  Range  for  Recognition  N=8  (Moser). 
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0  2000  4000  6000  8000  10000  12000  14000  16000  18000 

Range  in  meters 


N=  required  number  of  resels  on  critical  dimension  for  identification 
Figure  23  -  MRT  vs  Range  for  Identification  N=20  (Moser). 
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Temperature  difference  in 


MRT  vs  Range  for  Detection  N=6  (O'Neill) 


N=  required  number  of  resels  on  critical  dimension  for  detection 
Figure  24  -  MRT  vs  Range  for  Detection  N=6  (O’Neill). 
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0  0.5  1  1.5  2  2.5  3  3.5 

Range  in  meters  x  104 

N=  required  number  of  resels  on  critical  dimension  for  recognition 
Figure  25  -  MRT  vs  Range  for  Recognition  N=10  (O’Neill). 
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20 


MRT  vs  Range  for  Identification  N=32  (O'Neill) 


N=  required  number  of  resels  on  critical  dimension  for  identification 
Figure  26  -  MRT  vs  Range  for  Identification  N=32  (O’Neill). 
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V.  RESULTS 


After  the  procedures  described  in  the  previous  chapters,  the  sea  surface  radiance 
for  the  conditions  of  the  image  was  calculated  using  SEARAD  in  radiance  mode  for 
various  positions  of  the  FLIR  with  respect  to  the  target,  for  both  scenarios.  The  radiance 
values  and  the  MRT  values  were  tabulated  and  plotted  using  the  EXCEL  spreadsheet 
program.  The  plots  for  the  100  m  and  the  1000  m  sensor  elevations  are  displayed  in 
Figure  27  and  Figure  28  respectively. 

These  figures  show  the  Apparent  Target  to  Background  Temperature  Difference, 
ATapp,  with  no  polarizer  compared  with  that  for  horizontal  polarizer,  plotted  against  slant 
range  for  the  sensor  at  each  of  the  selected  fixed  elevations,  100  m  and  1000  m.  This 
represents  the  available  temperature  difference  at  the  receiver  as  a  function  of  range,  R, 
from  the  source.  The  Minimum  Resolvable  Temperature  Difference,  MRT,  is  shown  on 
the  same  scale  for  the  same  conditions.  The  intersection  of  the  available  AT  with  the 
MRT  defines  the  maximum  range  at  which  the  task  can  be  performed.  In  Figure  27  the 
MRT  for  identification  is  shown  for  sensor  at  100m  height.  The  range  at  which  this  MRT 
crosses  the  Apparent  Temperature  Difference  curve  defines  the  maximum  range  for 
identification.  Figure  28  shows  the  equivalent  information  for  sensor  at  1000m  height. 

For  the  100m  case  and  for  N=12  (six  line  pairs,  which  is  the  identification  case  for 
the  Johnson  criteria),  the  apparent  temperature  difference  for  the  horizontally  polarized 
case  crosses  the  MRT  line  at  a  range  of  13  km.  For  the  same  conditions  but  the 
unpolarized  case  the  two  lines  cross  at  a  range  of  about  28  km. 


83 


For  the  1000m  case  and  for  N=12  respectively,  the  apparent  temperature 
difference  for  the  horizontally  polarized  case  crosses  the  MRT  line  at  a  range  of  18.5  km. 

For  the  same  conditions  but  the  unpolarized  case  the  two  lines  meet  at  a  range  of 
about  29  km. 

The  inference  that  may  be  drawn  from  these  two  figures  is  that  the  use  of  the 
horizontal  polarizing  filter  decreases  the  maximum  range  for  identification,  although  the 
short-range  radiance  or  temperature  contrast  is  improved  by  the  filter.  This  loss  of 
identification  range  follows  from  the  fact  that  the  plotted  curves  of  ATapp  cross  at  a 
relatively  short  range,  and  beyond  this  “cross-over  range”  identification  range  is  reduced 
by  the  presence  of  the  filter.  This  conclusion  should  be  stated  with  caution  pending  a 
more  detailed  and  complete  study.  It  should  be  remembered  that  no  consideration  is  given 
to  reflected  sky  radiation  or  sun  glint  in  the  ship  signature,  with  or  without  polarization. 
The  possibility  of  error  in  the  magnitude  of  the  ship  signature  should  be  kept  in  mind, 
which  could  cause  error  in  the  apparent  temperature  difference,  and  hence  in  the  cross 
over  point. 

Two  anomalies  appear  in  the  curves  of  Figures  27  and  28,  namely; 

a)  From  the  plots  of  Figures  27  and  28  it  is  also  observed  that  at  some 
distance  beyond  the  crossover  range  the  temperature  difference  passes  through 
zero  and  goes  negative  for  the  horizontally  polarized  case 
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b)  In  Figure  28  both  the  unpolarized  and  the  polarized  signatures  exhibit  an 
apparent  upward  displacement  or  "bump"  in  the  range  around  6  to  10  km. 

At  first  glance  there  is  no  physical  explanation  to  this  behavior,  which  appears  to  be 
inherent  in  the  model  used. 

Since  at  greater  ranges  the  received  radiance  is  known  to  be  dominated  by  the 
path  radiance  which  is  common  to  both  source  and  background  signatures,  the  apparent 
temperatures  of  the  target  and  background  should  approach  each  other  asymptotically,  so 
that  the  apparent  temperature  difference  should  approach  zero.  The  unpolarized  apparent 
temperature  difference  shows  this  expected  behavior  in  both  scenarios,  but  in  both  cases 
the  horizontally  polarized  case  shows  this  anomaly.  This  suggests  that  the  model  used  is 
inadequate  in  some  way  to  describe  the  real  behavior.  Operationally,  the  anomaly  could 
be  made  to  disappear  by  adding  an  offset  of  about  1 1°  C  to  the  Horizontally  Polarized  AT 
curve  in  Figure  27,  and  a  similar  offset  to  Figure  28.  This  correction  is  displayed  in 
Figure  30,  where  the  temperature  is  forced  to  zero  at  36  km.  While  the  physical 
mechanism  of  such  an  offset  has  not  yet  been  identified,  several  hypotheses  have  been 
considered.  Additionally,  an  unexpectedly  low  value  of  unpolarized  temperature 
difference  occurs  at  very  short  ranges  in  the  1000m  scenario  only. 

The  negative  temperature  difference  might  perhaps  be  explained  by  considering 
the  different  spectral  distributions  of  the  target  and  the  background,  compared  with  the 
spectral  distribution  of  transmittance.  This  point  is  further  discussed  in  Chapter  VI, 
Discussion  and  Conclusions.  The  polarizers  for  the  wavelengths  of  interest  have 


85 


relatively  uniform  properties  over  the  band  used,  so  it  does  not  seem  that  they  should 
affect  adversely  the  spectral  distribution  of  the  target  and  the  background. 

The  increase  (bump)  in  the  apparent  temperature  difference  curve,  noted 
profoundly  in  the  unpolarized  case  for  sensor  at  1000m  at  ranges  up  to  8  km  is  of  interest. 
This  behavior  was  also  noted  in  Chih-Li  Yu  [Ref.  27].  It  may  have  a  connection  with  the 
modeling  of  the  waves  and  the  geometry  of  the  simulation  (including  Sun  direction),  but 
for  valid  conclusions  a  more  thorough  study  is  necessary.  This  is  also  considered  under 
discussion  in  the  next  chapter. 
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Difference  ii  °C 


FUR  at  100m  height,  N=12 


N=  required  number  of  resels  on  critical  dimension  for  identification 


Figure  27  -  MRT,  unpolarized  and  horizontally  polarized  apparent  temperature  difference  for 

sensor  at  100  meters. 


FUR  at  1000m  height,  N=12 


N=  required  number  of  resels  on  critical  dimension  for  identification 


Figure  28  -  MRT,  unpolarized  and  horizontally  polarized  apparent  temperature  difference  for 

sensor  at  1000  meters. 
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Sensor  at  1000m  heigth  Unpolarized 


1.07  1.13  1.91  2.5  2.99  3.32  3.98  4.6  5.46  6.73  8.79  10.4  12.7  16.4  17.2  17.7  18.3  18.9  19.6  20.3 

Slant  range  in  km 


Sea  Emission  Sky  Reflection  Sun  Glint 


Figure  29  -  Plot  of  full  range  (slant  ranges)  horizontally  polarized  values  from  SEARAD  for 
Sea  Emission,  Sky  Reflection  and  Sun  Glint  for  sensor  at  1000  meters  height. 
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FUR  at  100m  height,  N=12 


Range  in  Km 

»  MRT 1 0Om  — * —  Apparent  DT  Unp  -  -  Offset  App  DT  Hor  Pol  - Poly.  (Apparent  DT  Unp)  - Poly.  (Offset  App  DT  Hor  Pol) 


Figure  30  -  Plot  of  MRT,  unpolarized  and  horizontally  polarized  apparent  temperature 
difference  for  sensor  at  100  meters,  with  an  offset  of  1 1  degrees  Celsius  to  the  horizontally 

polarized  apparent  temperature  difference. 
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VI.  CONCLUSIONS  AND  DISCUSION 


From  the  plots  in  Figures  27  and  28  the  polarized  case  gives  a  significantly  bigger 
Apparent  Temperature  Difference  at  close  ranges  whereas  the  unpolarized  case  seems  to 
give  higher  Temperature  Differences  at  greater  ranges.  So  there  appears  to  be  a 
considerable  increase  in  the  Apparent  Temperature  Difference  by  using  a  horizontally 
polarized  filter  for  ranges  up  to  about  10  km.  From  the  outputs  of  this  work  there  are 
some  very  interesting  observations  which  require  interpretation. 

The  first  is  the  fact  that  the  Apparent  Temperature  Difference  for  the  polarized 
case  passes  through  zero  and  goes  negative  after  a  certain  range.  At  this  time  there  is  no 
physical  explanation  for  this  phenomenon.  The  temperature  differences  would  be 
expected  to  approach  zero  as  the  separate  radiances  become  dominated  by  the  path 
radiance  at  long  range.  A  possible  cause  of  this  anomaly  might  be  differences  in  the 
spectral  characteristics  of  the  target  and  background,  since  the  SEARAD  transmittances 
used  for  determination  of  apparent  target  temperature  are  averages  across  the  in-band 
wavelengths  of  interest  (8-12  /jm  in  this  case).  Thus  the  gray  body  distribution  that  was 
used  for  target  and  background  at  different  temperatures  might  have  caused  suppression 
(due  to  averaging)  of  strong  emission  or  absorption  sub-bands  within  the  8-12  /lm  band. 
With  a  spectral  analysis  of  the  target  and  background  perhaps  the  different  spectral 
characteristics  might  indicate  sub-bands  of  prominent  target  or  background  emission 
bands,  as  are  known  to  be  existent  in  the  3-5  jjm  band  already.  If  the  different  spectral 
characteristics  of  target  and  background  were  found  to  be  true,  it  could  lead  to  significant 
target  detection  advantages,  by  exploiting  the  spectral  characteristics  of  specific  targets. 
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In  the  scenario  studied  this  could  lead  to  preferential  propagation  of  either  source  or 
background  radiance  over  the  other. 

In  a  first  approach  to  this  issue,  a  plot  of  the  transmittance  by  sub-band  of  the  8- 
12/imband  was  made  using  SeaRad.  An  obvious  variation  in  the  transmittances  can  be 
seen  from  Figure  3 1 .  This  fact  has  an  important  effect  on  the  final  apparent  temperature 
difference,  since  the  poor  transmittance  in  one  part  of  the  sub-band  affects  negatively  the 
total  transmittance  of  the  whole  band.  The  hypothesis  is  presented  that  the  integral  over 
the  bandwidth  of  spectral  radiance  multiplied  by  spectral  transmittance  may  differ  from 
the  product  of  the  integrated  in-band  radiance  times  the  average  transmittance.  This  could 
possibly  also  help  provide  an  interpretation  of  the  related  issue,  that  the  two  Apparent 
Temperature  Difference  curves  for  the  horizontally  polarized  and  the  unpolarized  cases 
cross  each  other.  A  further  effort  to  evaluate  this  hypothesis  was  made  by  computing 
separately  the  product  of  radiance  and  transmittance  by  sub-band.  This  was  done 
separately  for  the  source  temperatures  of  the  sea  and  the  target,  with  the  results  shown  in 
Appendix  K.  The  conclusion  drawn  from  this  analysis  is  that  the  blackbody  spectral 
distributions  for  the  temperatures  of  target  and  background  differ  somewhat  in  magnitude 
but  negligibly  in  shape.  The  correctly  sub-band-integrated  transmittance  at  range  25-km 
showed  no  difference  to  four  significant  figures  between  the  two  source  temperatures. 
However  this  analysis  did  not  take  account  of  reflected  sky  radiance  in  the  ship  signature, 
while  this  is  a  significant  part  of  the  background  signature.  This  omission,  required  by 
the  level  of  target  modeling,  implies  an  under-estimate  of  the  apparent  temperature 
difference.  Target  size  and  shape  effects  are  not  included  in  this  model.  Representation 
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of  the  ship  target  by  its  hot  spot  temperature  is  also  somewhat  unrealistic,  since  the  sensor 
Detector  Angular  Subtense  of  one  (1)  mrad  gives  a  footprint  of  25  meters  at  25  km, 
which  greatly  exceeds  the  hot  spot  width.  The  conversion  of  source  temperature  to  and 
from  radiance  should  be  checked  further  for  consistency. 

The  "bump”  in  the  Apparent  Temperature  Difference  curve  in  Figures  27  and  28 
at  ranges  up  to  8  km  approximately  is  a  second  anomaly.  This  behavior  was  also  noted  in 
Chih-Li  Yu  [Ref.27].  In  these  ranges  the  zenith  angle  is  relatively  small,  approaching  the 
vertical.  In  Figure  29  the  separated  plots  of  the  full  range  (slant  ranges)  horizontally 
polarized  values  from  SEARAD  for  Sea  Emission,  Sky  Reflection  and  Sun  Glint  for 
sensor  at  1000  meters  height  reveal  that  the  Sky  Reflection  parameter  introduces  a  small 
bump  at  these  ranges.  It  is  also  known  from  Wolfe  and  Zissis  [Ref.26]  that  there  is  a 
variation  in  the  spectral  radiance  of  the  sky  at  different  angles  of  elevation.  These  facts 
along  with  the  wave  geometry  could  be  the  main  contributing  reason  for  this 
phenomenon.  These  are  the  first  thoughts  about  this  and  more  study  is  needed  to  derive 
more  complete  conclusions. 

Conclusions  as  to  the  relative  values  of  identification  range  for  polarized  and 
unpolarized  cases  should  be  treated  as  tentative  pending  further  analysis  and 
improvement  of  the  model. 
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Transmittanc 


8-12  microns  for  FLIR  at  100m  height  and  24.5  Km  range 


Wavelength  in  microns 

Figure  3 1  -  Transmittances  in  relation  to  the  sub-band  wavelength. 
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VII.  RECOMMENDATIONS  FOR  FUTURE  WORK 


One  observation  of  interest  in  this  work  is  the  negative  apparent  temperature 
difference  of  the  horizontally  polarized  case  after  a  certain  range.  Due  to  the  lack  of  any 
physical  explanation  of  the  phenomenon  at  this  time  some  new  ideas  about  the  target  and 
background  characteristics  were  considered.  One  such  explanation  deriving  from  the  3- 
5  jjm  band  relates  with  the  different  spectral  properties  of  target  and  background. 

The  target  and  background  behave  in  a  different  way  spectrally.  In  this 
experiment  SeaRad  calculated  average  values  for  the  wavelengths  of  interest  (8- 
12.05 fim).  In  this  way  a  wavelength  region  of  low  transmittance  of  the  target  could 
effect  adversely  the  infrared  characteristics  of  the  target.  Analyzing  that  in  the  spectral 
domain  and  examining  the  properties  of  target  and  background  spectrally  could  provide 
characteristic  patterns  that  would  eventually  improve  discrimination  of  the  target  from 
the  background. 

It  is  realized  that  the  standard  formulation  of  the  MRT  is  based  on  a  300°  K 
background.  Since  the  background  temperature  varies  and  is  only  sometimes  that  the 
specific  300°  K  matches  for  the  real  conditions,  a  further  study  of  the  MRT  formulation  is 
needed  in  order  to  examine  the  accuracy  of  the  calculations  in  other  cases  than  the  300°  K 
background. 

The  algorithm  used  for  conversion  of  source  temperature  to  radiance,  and  its 
inversion  to  return  to  apparent  temperature  should  be  evaluated  and  validated, 
particularly  as  applied  to  partially  polarized  radiation. 
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APPENDIX  A 


EXPERIMENT  LOCATION 

Description : 

Plots  of  the  GPS  information  of  the  ship  trajectory  during  the  experiment 
in  April  9  and  April  10,  1996.  Image  set  recording  points  are  marked  as  diamonds  in  this 
figure. 

Map  showing  the  planned  ship  positioning  and  desired  ship  heading  in  each 
station  location. 

Conditions : 

The  AGA-780  Thermovision  was  located  at  the  Building  15  of  the 
NCCOSC-NRaD  in  Point  Loma.  Its  geographic  coordinates  were  N  32  39  35.8,  W  117 
14  22.2.  It  was  mounted  on  a  table  projecting  from  the  window,  and  directed  using  a 
remote  pan-tilt  head.  The  ship  maneuvered  in  reference  to  the  camera  position  according 
to  five  planned  magnetic  bearings  (camera  to  ship)  ranging  from  135  to  200  degrees.  For 
each  bearing  there  are  5  station  locations  at  ranges  from  0.5  to  3.0  nm.  In  each  station 
location  the  ship  is  supposed  to  turn  to  the  indicated  magnetic  headings  and  stay  still 
during  the  time  needed  for  the  image  recording  process.  This  process  was  modified  by 
ship-handling  necessities. 
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Notes : 


Both  ship  GPS  position  and  camera  location  are  known,  so  bearing  is  computed. 
Thus,  ship  heading  information  is  used  with  the  bearing  information  to  give  the  aspect 
angle  of  the  target. 

Range  to  the  shore  station  and  ship  heading  were  reported  from  the  ship  at 
irregular  intervals  and  recorded  at  the  sensor.  Magnetic  bearings  to  the  ship  were 
recorded  on  land  in  the  sensor  log. 
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Experiment  area  (San  Diego,  California).  The  courses  on  the  map  are  the  courses  that 
R/V  POINT  SUR  was  to  follow  during  the  experiment. 
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Ship’s  trajectory  on  9  April  1996. 
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Ship’s  trajectory  on  10  April  1996. 
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APPENDIX  B 


POLARIZER  FILTERS 

Description : 

The  two  following  figures  show  the  performance  curves  for  the  KRS-5  substrate, 
9.5mm  polarizers  used  in  the  camera  setting  positions  4  and  5  respectively. 

Conditions : 

The  data  were  provided  by  the  manufacturer  of  the  polarizers,  Graseby-Specac. 
The  measurements  were  taken  with  a  spectrometer  using  the  actual  filter  shipped.  Where 
two  filters  were  required  for  the  crossed-grid  measurement  (trace  3),  another  polarizer 
nearly  identical  to  the  one  shipped  is  used. 

Notes : 

Trace  1:  Ev  i.e.  optical  power  transmitted  by  the  polarizer  with  its  grid  lines 
vertical. 

Trace  2:  Eh  i.e.  optical  power  transmitted  by  the  polarizer  with  its  grid  lines 
horizontal. 

Trace  3:  Ei  i.e.  optical  power  transmitted  by  a  pair  of  polarizers  with  grids  crossed 
Vertical  axis  :  Wavenumber/wavelength 
Horizontal  axis:  Transmittance 
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The  principal  transmittances  Ki  and  K2  of  a  polarizer  are  defined  as  the  transmittances  of 
the  polarizer  for  radiation  whose  electric  field  vector  is  (respectively)  perpendicular  to 
and  parallel  to  the  grid  lines.  The  full  equations  for  Ki  and  K2  are: 


K,= 


2E0 


",  4£,£0 

1/2 

i 

CN 

feT 

+ 

1 

H 

- 1 

K,  =  Eh*E’ 


2  En 


1- 


1- 


4£,£0 

{Eh+Evf 


1/2 


In  these  equations,  Eo  is  the  total  power  transmitted  through  the  measurement  system 
with  no  polarizers  in  the  beam.  For  the  attached  traces,  the  spectrometer  output  was 
normalized  such  that  E0=  1.00  (i.e.  100%)  for  all  wavelengths.  As  can  be  seen  from  trace 
3,  Ei  is  small  over  the  operating  wavelength  range  of  the  polarizer.  Consequently,  the 
above  equation  may  be  simplified  (for  wavelengths  within  the  operating  range)  to  the 
following  approximate  expressions: 


_  Eh+  Ev 

E, 


K  2  = 


gi 

Eh  +  Ev 
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So  for  the  specific  case  the  bandwidth  of  interest  is  from  8  to  12.5  jum .  From  the  curves 
Eh=0. 17  and  Ev=0.53.  Eq=1,  thus 


.  _£,+£„  _  0.17  +  0.53 
'  ”  E0  1 


The  polarizing  filter  transmittance  for  non  polarized  waves  in  the  8  to  12.5  jum ,  is  0.7. 
For  this  region  the  average  value  from  the  Trace  2  (Eh  i.e.  optical  power  transmitted  by 
the  polarizer  with  its  grid  lines  horizontal)  is  0.5.  So  these  two  values  will  multiply  the 
background  and  the  ship  plus  path  radiance  respectively  to  provide  the  energy  reaching 
the  sensor  after  passing  through  the  polarizing  filter. 
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Research  Vessel  Point  Sur 


Location  of  the  thermistors  onboard  R/V  Point  Sur. 
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APPENDIX  D 


| 


SEARAD 

This  is  an  example  of  the  tape5-input  files  used  in  SeaRad.  The  atmospheric  model 


was  obtained  from  the  meteorological  data  from  the  experiment. 


7  3 

2  1 

2  2 

2  2 

3  1 

0  7 

0  0 

0.000 

34  0 

0 

METEOR  DATA 

0.003 

1018.000 

14.300 

75.000 

0.009 

1017.700 

15.000 

63.000 

0.017 

1016.700 

14.800 

63.000 

0.034 

1014.700 

14.300 

65.000 

0.046 

1013.200 

14.000 

66.000 

0.056 

1012.000 

13.900 

67.000 

0.069 

1010.500 

13.800 

68.000 

0.075 

1009.800 

13.700 

68.000 

0.087 

1008.300 

13.600 

69.000 

0.100 

1006.800 

13.500 

69.000 

0.112 

1005.300 

13.400 

69.000 

0.130 

1003.100 

13.200 

70.000 

0.155 

1000.200 

13.000 

71.000 

0.163 

999.200 

12.800 

72.000 

0.183 

996.800 

12.700 

72.000 

0.216 

992.900 

12.400 

73.000 

0.243 

989.800 

12.100 

74.000 

0.267 

986.900 

12.000 

75.000 

0.304 

982.500 

11.800 

76.000 

0.332 

979.200 

11.700 

76.000 

0.381 

973.500 

11.500 

77.000 

0.418 

969.200 

11.200 

78.000 

0.497 

960.000 

11.200 

77.000 

0.525 

956.700 

11.300 

77.000 

0.588 

949.500 

11.100 

77.000 

0.613 

946.700 

11.300 

74.000 

0.655 

941.900 

12.100 

56.000 

0.696 

937.300 

14.300 

24.000 

0.738 

932.700 

15.800 

12.000 

0.800 

925.900 

15.500 

12.000 

0.847 

920.800 

15.100 

12.000 

0.892 

915.800 

15.300 

9.000 

0.956 

908.900 

15.500 

8.000 

1.004 

903.800 

15.700 

7.000 

0.100 

0.000 

90.310 

00.000 

1  2 

99  0 

32.665 

117.242 

0.000 

0.000 

833 

1250 

10 

5 

o 


22110  289.95  -5 

7.000  4.481  0.000  0.000 

ABH 

ABH 

ABH 

ABH 

ABH 

ABH 

ABH 

ABH 

ABH 

ABH 

ABH 

ABH 

ABH 

ABH 

ABH 

ABH 

ABH 

ABH 

ABH 

ABH 

ABH 

ABH 

ABH 

ABH 

ABH 

ABH 

ABH 

ABH 

ABH 

ABH 

ABH 

ABH 

ABH 

ABH 

0.000  0.00  0  85.077  T 

20.238  198.923  0.000  0.000 

0 
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This  is  an  example  of  a  SeaRad  out  file  of  the  previous  tape5  file. 

*****  SEARAD,  A  MODIFICATION  OF  LOWTRAN7  ***** 

DATE:  03/02/1999 
TIME:  20:43:37.21 


THERMAL  PLUS  SOLAR  RADIANCE  MODE 
MULTIPLE  SCATTERING  USED 


MARINE  AEROSOL  MODEL  USED 


WIND  SPEED 
WIND  SPEED  = 

RELATIVE  HUMIDITY  = 
AIRMASS  CHARACTER  =  ' 
VISIBILITY 


7.00  M/SEC 

4.48  M/SEC,  24  HR  AVERAGE 
74.94  PERCENT 
7.0 

.00  KM 


SLANT  PATH  TO  SPACE 


HI 

= 

.100 

KM 

HMIN 

= 

.000 

KM 

ANGLE 

= 

90.310 

DEG 

FREQUENCY  RANGE 


IV1 

830 

CM-l 

{  12.05 

MICROMETERS) 

IV2  = 

1250 

CM-1 

(  8.00 

MICROMETERS) 

IDV  = 

10 

CM-l 

IFWHM 

5 

CM-1 

IFILTER  = 

0 

SUMMARY  OF  THE  GEOMETRY  CALCULATION 


HI 

= 

.100 

KM 

H2 

= 

.000 

KM 

ANGLE 

= 

90.310 

DEG 

RANGE 

= 

29.376 

KM 

BETA 

= 

.264 

DEG 

PHI 

= 

89.901 

DEG 

HMIN 

= 

.000 

KM 

BENDING 

= 

.053 

DEG 

LEN 

= 

0 

SEA  AT  289.95  K  REPLACES  BLACK  BODY  BOUNDARY 

UPWIND  =  85.077  DEG  EAST  OF  LINE  OF  SIGHT 


SUMMARY  OF  OBSERVATION  GEOMETRY 
BETA 

PATH  AZIMUTH 
WIND  AZIMUTH 

RECEIVER  LATITUDE 
RECEIVER  LONGITUDE 
FOOTPRINT  LATITUDE 
FOOTPRINT  LONGITUDE 
SUBSOLAR  LATITUDE 


=  .26417  DEG 

198.923  DEG  EAST  OF  NORTH 
284.000  DEG  EAST  OF  NORTH 

32.665  NORTH  OF  EQUATOR 
117.242  WEST  OF  GREENWICH 
32.415  NORTH  OF  EQUATOR 
117.343  WEST  OF  GREENWICH 
7.057  DEG  NORTH  OF  EQUATOR 
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SUBSOLAR  LONGITUDE 


122.989  DEG  WEST  OF  GREENWICH 


VALUES  SEEN  FROM  FOOTPRINT 


RECEIVER  ZENITH  ANGLE 

= 

89.901 

DEG 

RECEIVER  AZIMUTH 

= 

265.132 

DEG  WEST  OF 

UP  WIND 

SOLAR  ZENITH  ANGLE 

= 

25.897 

DEG 

SOLAR  AZIMUTH 

= 

91.084 

DEG  WEST  OF 

UP  WIND 

SOLAR  SPECULAR  TILT 

= 

32.192 

DEG  (  4.89 

SIGMA,  PROB= 

5 . 282E-05 ) 

ZERO 

RANGE  UNPOLARIZED  VALUES 

SEA  EMISSION 

= 

20.01836 

W  M-2  SR-1 

(AV.  EMISS. 

.6071) 

SKY  REFLECTION 

= 

4.61739 

W  M-2  SR-1 

SUN  GLINT 

= 

.00002  W 

M-2  SR-1 

TOTAL  RADIANCE 

= 

24.63576 

W  M-2  SR-1 

BLACK  BODY  TEMP. 

= 

1.1 

C 

FULL 

RANGE  UNPOLARIZED  VALUES 

SEA  EMISSION 

= 

1.02095 

W  M-2  SR-1 

SKY  REFLECTION 

= 

.21319 

W  M-2  SR-1 

SUN  GLINT 

= 

.00000  W 

M-2  SR-1 

PATH  TO  FOOTPRINT 

= 

29.24329 

W  M-2  SR-1 

(AV.  TRANS. 

.0514) 

TOTAL  RADIANCE 

- 

30.47743 

W  M-2  SR-1 

BLACK  BODY  TEMP. 

= 

12.5 

C 

HORIZONTAL  VERTICAL 

(H-V)  /  (H+V) 

(W  M- 

2  SR-1)  (W  M-2  SR-1) 

(%) 

ZERO 

RANGE  POLARIZED  VALUES 

SEA  EMISSION 

8.88451 

11.13385 

-11.2 

SKY  REFLECTION 

2.63338 

1.98401 

14.1 

SUN  GLINT 

.00002 

.00000 

92.0 

TOTAL  RADIANCE 

11.51791 

13.11785 

-6.5 

BLACK  BODY  TEMP.  <C) 

- 

33.5 

-28.2 

FULL 

RANGE  POLARIZED  VALUES 

SEA  EMISSION 

.45213 

.56882 

-11.4 

SKY  REFLECTION 

.12153 

.09166 

14.0 

SUN  GLINT 

.00000 

.00000 

.0 

PATH  TO  FOOTPRINT 

14, 

.62165 

14.62165 

.0 

TOTAL  RADIANCE 

15, 

.19530 

15.28213 

-.3 

BLACK  BODY  TEMP.  (C) 

-21.9 

-21.6 
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APPENDIX  E 


METEOROLOGICAL  DATA 

1.  METOC  FILES 

Storage  format  for  METOCOl  (DASC)  data 

The  data  is  stored  in  a  set  of  files  with  the  names  "eopace.*"  and  "30sec.*".  The 
extension  for  the  file  name  is  the  Julian  data  the  file  was  started  on.  A  new  pair  of  files  is 
started  every  day  at  00:00  GMT.  The  time  data  in  the  files  are  in  GMT.  The  eopace  files 
have  the  data  averaged  over  10  minutes.  The  30sec  files  have  the  GPS  data  averaged  over 
30  seconds  and  the  last  set  of  data  from  the  Campbell.  In  the  30sec  file,  the  Campbell 
data  are  not  averaged.  The  data  format  for  both  files  is  the  same.  Each  file  has  a  header 
that  tells  when  the  file  was  started.  Each  file  has  a  footer  that  tells  when  the  file  was 
finished.  Some  files  may  also  have  additional  start  and  stop  times  interspersed  with  data, 
depending  on  whether  data  logging  was  turned  off  or  the  program  was  stopped  before  the 
end  of  the  day.  The  data  are  arranged  in  fields  separated  by  spaces  across  a  line  of  text.  A 
carriage  return  character  and  a  line  feed  character  are  at  the  end  of  each  line.  Each  new 
reading  is  on  a  new  line.  The  fields  are  in  the  order  below.  The  numbers  before  the  fields 
are  not  in  the  data.  They  serve  to  make  identifying  fields  easier. 

01  Number  of  readings  in  this  average,  nn 

02  Year,  yyyy 

03  Julian  date,  ddd 

04  Hours  and  minutes,  hhmm 
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05  Seconds,  ss 

06  Relative  wind  speed,  m/s,  mmm 

07  Relative  wind  direction,  degrees,  ddd 

08  Ship  speed,  knots,  at  present  always  0 

09  Ship  direction,  degrees,  at  present  always  090 

10  True  wind  speed,  m/s,  mmm 

1 1  True  wind  direction,  degrees,  ddd 

12  T  air,  degrees  C,  tt.t 

13  RH,  per  cent,  pp 

14  Pressure,  millibars,  pppp.p 

1 5  Sea  surface  temperature,  degrees  C,  tt.t 

16  GPS  time,  hhmmss 

17  GPS  latitude,  ddmm.m 

1 8  GPS  North/South  indicator,  either  ’N’  or  ’S’ 

19  GPS  longitude,  dddmm.m 

20  GPS  East/West  indicator,  either  E’  or  W’ 

21  GPS  speed  over  ground,  knots,  kk.k 

22  GPS  course  over  ground,  degrees,  ddd. 

23  GPS  antenna  height,  meters,  mm.m 

An  example  of  the  storage  format  is  presented  in  the  following  pages. 
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Storage  format  for  METOC02  data 

The  METOC02  system  had  two  additionsl  temperature  sensors  and  two  additional 
RH  sensors.  The  data  are  arranged  slightly  diffemtly  than  the  METOCOl  system.  The 
data  is  stored  in  a  set  of  files  with  the  names  "eopace.*"  and  "30sec.*".  The  extension  for 
the  file  name  is  the  Julian  data  the  file  was  started  on.  A  new  pair  of  files  is  started  every 
day  at  00:00  GMT.  The  time  data  in  the  files  are  in  GMT.  The  eopace  files  have  the  data 
averaged  over  10  minutes. 

The  30sec  file  have  the  GPS  data  averaged  over  30  seconds  and  the  last  set  of  data 
from  the  Campbell.  In  the  30sec  file,  the  Campbell  data  are  not  averaged.  The  data  format 
for  both  files  is  the  same.  Each  file  has  a  header  that  tells  when  the  file  was  started.  Each 
file  has  a  footer  that  tells  when  the  file  was  finished.  Some  files  may  also  have  additional 
start  and  stop  times  interspersed  with  data,  depending  on  whether  data  logging  was  turned 
off  or  the  program  was  stopped  before  the  end  of  the  day.  The  data  are  arranged  in  fields 
separated  by  spaces  across  a  line  of  text.  A  carriage  return  character  and  a  line  feed 
character  are  at  the  end  of  each  line.  Each  new  reading  is  on  a  new  line.  The  fields  are  in 
the  order  below.  The  numbers  before  the  fields  are  not  in  the  data.  They  serve  to  make 
identifying  fields  easier. 

01  Number  of  readings  in  this  average,  nn 

02  Year,  yyyy 

03  Julian  date,  ddd 

04  Hours  and  minutes,  hhmm 

05  Seconds,  ss 
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06 

07 

08 

09 

10 

11 

12 

13 

14 

15 

16 

17 

18 

19 

20 
21 
22 

23 

24 

25 

26 
27 


Relative  wind  speed,  m/s,  mmm 

Relative  wind  direction,  degrees,  ddd 

Ship  speed,  knots,  at  present  always  0 

Ship  direction,  degrees,  at  present  always  090 

True  wind  speed,  m/s,  mmm 

True  wind  direction,  degrees,  ddd 

T  air,  degrees  C,  tt.t 

RH,  per  cent,  pp 

Pressure,  millibars,  pppp.p 

Sea  surface  temperature,  degrees  C,  tt.t 

T  air  A,  degrees  C,  tt.t 

RH  A,  per  cent,  pp 

T  air  B,  degrees  C,  tt.t 

RH  B,  per  cent,  pp 

Delimiter  between  Campbell  data  and  GPS  data,  "***" 

GPS  time,  hhmmss 

GPS  latitude,  ddmm.m 

GPS  North/South  indicator,  either  ’N’or  ’S’ 

GPS  longitude,  dddmm.m 

GPS  East/West  indicator,  either  E’  or  W’ 

GPS  speed  over  ground,  knots,  kk.k 
GPS  course  over  ground,  degrees,  ddd. 
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28  GPS  antenna  height,  meters,  mm.m 

An  example  of  the  storage  format  is  presented  in  the  following  page. 
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METOC  1  FORMAT 


Data  logging  started  at  Tue  Apr  09  00:00:00  1996 

30  1996  99  2355  44  2  156  0  090  6  294  15.4  81  1016.2 

21.2  000024  3311.2  N  11728.2  W  6.3  124.3  56.8 

30  1996  99  2355  44  2  156  0  090  6  294  15.4  81  1016.2 

21.2  000054  3311.2  N  11728.2  W  6.3  123.9  52.5 

30  1996  99  2355  44  2  156  0  090  6  294  15.4  81  1016.2 

21.2  000124  3311.1  N  11728.1  W  6.3  123.5  42.0 

30  1996  99  2355  44  2  156  0  090  6  294  15.4  81  1016.2 

21.2  000154  3311.1  N  11728.0  W  6.3  123.0  26.9 

30  1996  99  2355  44  2  156  0  090  6  294  15.4  81  1016.2 

21.2  000224  3311.1  N  11728.0  W  6.4  121.6  9.9 

30  1996  99  2355  44  2  156  0  090  6  294  15.4  81  1016.2 

21.2  000254  3311.1  N  11727.9  W  6.3  120.5  -4.2 

30  1996  99  2355  44  2  156  0  090  6  294  15.4  81  1016.2 

21.2  000324  3311.0  N  11727.9  W  5.9  129.6  -10.7 

30  1996  99  2355  44  2  156  0  090  6  294  15.4  81  1016.2 

21.2  000354  3311.0  N  11727.9  W  4.1  185.6  -6.7 

30  1996  99  2355  44  2  156  0  090  6  294  15.4  81  1016.2 

21.2  000424  3311.0  N  11727.9  W  4.1  233.5  9.0 


METOC  2  FORMAT 


Data  logging  started  at  Tue  Apr  09  00:00:00  1996 

30  1996  99  2352  32  2  138  0  090  5 

288 

16.5 

75  1016.9 

16.8 

0.0  0.0  0.0  0.0***  000035 

3311.2  N  11728.2  W 

6.3 

124.3 

46.3 

30 

1996  99  2352  32  2 

138  0  090  5 

288 

16.5 

75  1016.9 

16.8 

0.0  0.0  0.0  0.0***  000105 

3311.1  N  11728.1  W 

6.3 

124.2 

40.7 

30 

1996  99  2352  32  2 

138  0  090  5 

288 

16.5 

75  1016.9 

16.8 

0.0  0.0  0.0  0.0***  000135 

3311.1  N  11728.1  W 

6.2 

123.6 

29.4 

29 

1996  99  2352  32  2 

138  0  090  5 

288 

16.5 

75  1016.9 

16.8 

0.0  0.0  0.0  0.0***  000205 

3311.1  N  11728.0  W 

6.4 

122.3 

14.0 

30 

1996  100  0002  8  2 

149  0  090  5 

292 

16.2 

76  1016.9 

16.6 

0.0  0.0  0.0  0.0***  000235 

3311.1  N  11728.0  W 

6.3 

121.4 

-2.1 

30 

1996  100  0002  8  2 

149  0  090  5 

292 

16.2 

76  1016.9 

16.6 

0.0  0.0  0.0  0.0***  000305 

3311.0  N  11727.9  W 

6.2 

120.7 

-13.8 

30 

1996  100  0002  8  2 

149  0  090  5 

292 

16.2 

76  1016.9 

16.6 

0.0  0.0  0.0  0.0***  000335 

3311.0  N  11727.9  W 

4.6 

159.8 

-17.5 

30 

1996  100  0002  8  2 

149  0  090  5 

292 

16.2 

76  1016.9 

16.6 

0.0  0.0  0.0  0.0***  000405 

3311.0  N  11727.9  W 

3.7 

216.5 

-10.2 

30 

1996  100  0002  8  2 

149  0  090  5 

292 

16.2 

76  1016.9 

16.6 

0.0  0.0  0.0  0.0***  000435 

3311.0  N  11727.9  W 

3.8 

259.1 

8.3 
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2.  RADIOSONDE  DATA 

The  following  data  are  the  radiosonde  data  collected  during  the  experiment.  The  data 
were  provided  by  the  Meteorological  department  of  the  Naval  Postgraduate  School. 

Sounding  program  REV  7.62  using  Omega 

Ship  :  POINT  SUR 

Location  :  32.95  N  117.44  W  3m 

Rejected  Sigma  stations:  a,b,c 

Phase  fitting  length  is  250  s  from  0  min  to  120  min 


Sounding  :  44 
RS-number:  542520109 


No  PTU  editing 
No  wind  editing 

Started  at:  9  APR  96  13:49  GMT 


Time  AscRate  Hgt/MSL  Pressure  Temp  RH  Dewp  Dir  Speed  WndStat 


min 

s 

m/s 

m 

hPa 

degC 

% 

degC 

deg 

m/s 

0 

0 

0.0 

3 

1017.9 

14.0 

74 

9.5 

279 

3.3 

0 

2 

4.0 

11 

1016.9 

15.3 

62 

8.1 

m 

nn 

--CDE-GH- 

0 

4 

4.7 

22 

1015.7 

15.2 

61 

7.8 

m 

nn 

— CDE-GH- 

0 

6 

6.2 

40 

1013.5 

14.8 

64 

8.1 

m 

nn 

--CDE-GH- 

0 

8 

4.9 

42 

1013.2 

14.4 

64 

7.7 

m 

n  /  / 

—CDE-GH- 

0 

10 

5.1 

54 

1011.7 

14.1 

65 

7.7 

328 

1.4 

--CDE-GH- 

0 

12 

5.2 

65 

1010.5 

13.8 

66 

7.6 

m 

nn 

—CDE-GH- 

0 

14 

5.0 

73 

1009.5 

13.7 

67 

7.7 

m 

/ 1  /  / 

—CDE-GH- 

0 

16 

5.1 

85 

1008.0 

13.6 

69 

8.1 

m 

//// 

— CDE-GH- 

0 

18 

4.8 

89 

1007.6 

13.5 

69 

8.0 

m 

//// 

--CDE-GH- 

0 

20 

4.8 

100 

1006.3 

13.4 

70 

8.1 

343 

0.9 

--CDE-GH- 

0 

22 

5.0 

114 

1004.6 

13.3 

70 

8.0 

m 

nn 

--CDE-GH- 

0 

24 

4.8 

118 

1004.1 

13.2 

70 

7.9 

m 

/  n  / 

— CDE-GH- 

0 

26 

5.1 

136 

1001.9 

13.1 

71 

8.0 

m 

//// 

— CDEFGH- 

0 

28 

4.7 

13  6 

1001.9 

13.0 

71 

7.9 

m 

//// 

--CDEFGH- 

0 

30 

4.7 

143 

1001.2 

13.0 

71 

7.9 

347 

1.2 

--CDEFGH- 

0 

32 

4.9 

159 

999.2 

12.8 

72 

7.9 

III 

till 

— CDEFGH- 

0 

34 

4.7 

163 

998.7 

12.8 

73 

8.1 

III 

nn 

--CDEFGH- 

0 

36 

4.7 

181 

996.5 

12.6 

73 

7.9 

III 

nn 

— CDEFGH- 

0 

38 

4.8 

187 

995.8 

12.6 

74 

8.1 

III 

n  n 

--CDEFGH- 

0 

40 

4.8 

198 

994.6 

12.5 

74 

8.0 

7 

0.8 

--CDEFGH- 

0 

42 

4.8 

210 

993.1 

12.3 

74 

7.9 

III 

//// 

--CDEFGH- 

0 

44 

4.8 

216 

992.4 

12.3 

74 

7.9 

III 

//// 

--CDEFGH- 

0 

46 

4.8 

230 

990.7 

12.1 

75 

7.9 

III 

//// 

--CDEFGH- 

0 

48 

4.8 

234 

990.2 

12.0 

76 

8.0 

III 

//// 

— CDEFGH- 

0 

50 

4.8 

245 

989.0 

12.0 

76 

8.0 

38 

0.5 

— CDEFGH- 

123 


0 

52 

00 

259 

987.3 

11.8 

77 

o 

00 

m 

/III 

--CDEFGH- 

0 

54 

5.0 

269 

986.1 

11.7 

77 

7.9 

m 

llll 

--CDEFGH- 

0 

56 

5.0 

287 

984.0 

11.6 

78 

7.9 

m 

llll 

--CDEFGH- 

0 

58 

5.2 

291 

983.5 

11.5 

78 

7.8 

m 

llll 

--CDEFGH- 

1 

0 

5.2 

300 

982.5 

11.4 

78 

7.8 

32 

0.5 

--CDEFGH- 

1 

2 

5.1 

312 

981.1 

11.3 

79 

7.8 

III 

llll 

--CDEFGH- 

1 

4 

5.2 

318 

980.4 

11.2 

79 

7.7 

III 

llll 

--CDEFGH- 

1 

6 

4.8 

324 

979.7 

11.1 

79 

7.6 

III 

llll 

--CDEFGH- 

1 

8 

4.8 

330 

978.9 

11.1 

80 

7.8 

III 

llll 

--CDEFGH- 

1 

10 

4.7 

340 

977.7 

11.0 

80 

7.7 

0 

0.0 

--CDEFGH- 

1 

12 

4.7 

350 

976.6 

10.9 

81 

7.8 

III 

llll 

--CDEFGH- 

1 

14 

4.8 

359 

975.6 

10.8 

81 

7.7 

III 

llll 

--CDEFGH- 

1 

16 

4.6 

369 

974.4 

10.8 

81 

7.7 

III 

llll 

--CDEFGH- 

1 

18 

4.8 

379 

973.2 

10.7 

82 

7.8 

III 

llll 

--CDEFGH- 

1 

20 

4.7 

385 

972.5 

10.6 

82 

7.7 

0 

0.0 

--CDEFGH- 

1 

22 

4.6 

397 

971.1 

10.3 

83 

7.6 

III 

llll 

--CDEFGH- 

1 

24 

4.4 

401 

970.6 

10.4 

83 

7.7 

III 

llll 

--CDEFGH- 

1 

26 

4.1 

411 

969.4 

10.3 

84 

7.8 

III 

llll 

--CDEFGH- 

1 

28 

4.3 

419 

968.5 

10.2 

84 

7.7 

III 

llll 

--CDEFGH- 

1 

30 

4.3 

429 

967.3 

10.1 

85 

7.7 

0 

0.0 

--CDEFGH- 

1 

32 

4.3 

440 

966.1 

10.0 

85 

7.6 

III 

llll 

--CDEFGH- 

1 

34 

4.3 

446 

965.4 

10.0 

86 

7.8 

III 

llll 

--CDEFGH- 

1 

36 

4.4 

456 

964.2 

9.8 

87 

7.8 

III 

llll 

--CDEFGH- 

1 

38 

4.3 

460 

963.8 

9.8 

87 

7.8 

III 

llll 

--CDEFGH- 

1 

40 

4.2 

466 

963.1 

9.7 

88 

7.9 

0 

0.0 

--CDEFGH- 

1 

42 

4.3 

478 

961.7 

9.6 

88 

7.8 

III 

llll 

--CDEFGH- 

1 

44 

4.2 

486 

960.7 

9.6 

89 

7.9 

III 

llll 

--CDEFGH- 

1 

46 

4.1 

492 

960.0 

9.5 

89 

7.8 

III 

llll 

--CDEFGH- 

1 

48 

4.3 

508 

958.1 

9.4 

90 

7.9 

III 

llll 

--CDEFGH- 

1 

50 

4.3 

514 

957.4 

9.3 

90 

7.8 

178 

0.6 

— CDEFGH- 

1 

52 

4.5 

533 

955.3 

9.1 

91 

7.8 

III 

llll 

--CDEFGH- 

1 

54 

4.7 

541 

954.4 

9.1 

91 

7.8 

III 

/III 

--CDEFGH- 

1 

56 

4.5 

545 

953.9 

9.0 

91 

7.7 

III 

llll 

--CDEFGH- 

1 

58 

4.7 

559 

952.3 

8.9 

92 

7.7 

III 

llll 

--CDEFGH- 

2 

0 

4.5 

563 

951.8 

8.8 

92 

7.6 

182 

0.8 

--CDEFGH- 

2 

2 

4.4 

573 

950.7 

8.7 

93 

7.7 

III 

llll 

--CDEFGH- 

2 

4 

4.5 

581 

949.8 

8.6 

93 

7.6 

III 

llll 

--CDEFGH- 

2 

6 

4.3 

585 

949.3 

8.6 

93 

7.6 

III 

llll 

--CDEFGH- 

2 

8 

4.6 

597 

947.9 

8.5 

94 

7.6 

III 

llll 

— CDEFGH- 

2 

10 

4.6 

605 

947.0 

8.4 

94 

7.5 

181 

1.1 

--CDEFGH- 

2 

12 

4.9 

625 

944.7 

8.3 

95 

7.6 

III 

llll 

--CDEFGH- 

2 

14 

5.0 

637 

943.3 

8.2 

95 

7.5 

III 

llll 

--CDEFGH- 

2 

16 

5.0 

641 

942.8 

8.1 

96 

7.5 

III 

till 

--CDEFGH- 

2 

18 

5.1 

661 

940.5 

8.0 

96 

7.4 

III 

llll 

--CDEFGH- 

2 

20 

5.2 

671 

939.4 

8.0 

96 

7.4 

171 

1.7 

--CDEFGH- 

2 

22 

4.9 

681 

938.2 

7.9 

96 

7.3 

III 

llll 

--CDEFGH- 

2 

24 

5.0 

690 

937.3 

7.8 

96 

7.2 

III 

llll 

--CDEFGH- 

2 

26 

5.2 

700 

936.2 

7.7 

96 

7.1 

III 

llll 

--CDEFGH- 

2 

28 

5.2 

714 

934.6 

7.6 

96 

7.0 

III 

llll 

--CDEFGH- 

2 

30 

5.2 

720 

933.9 

7.5 

96 

6.9 

175 

2.0 

--CDEFGH- 

2 

32 

5.3 

732 

932.5 

7.4 

96 

6.8 

III 

llll 

--CDEFGH- 

2 

34 

5.2 

738 

931.8 

7.4 

96 

6.8 

III 

llll 

--CDEFGH- 

2 

36 

5.3 

744 

931.2 

7.4 

91 

6.1 

III 

llll 

--CDEFGH- 

2 

38 

5.4 

760 

929.3 

7.7 

88 

5.9 

III 

llll 

--CDEFGH- 

124 


2 

40 

5.4 

768 

928.4 

7.9 

82 

5.1 

176 

2.5 

--CDEFGH- 

2 

42 

5.3 

784 

926.6 

10.1 

42 

-2.2 

III 

llll 

--CDEFGH- 

2 

44 

5.1 

790 

925.9 

11.6 

28 

-6.2 

III 

llll 

--CDEFGH- 

2 

46 

5.3 

800 

924.8 

12.6 

20 

-9.7 

III 

nn 

--CDEFGH- 

2 

48 

5.0 

810 

923.7 

14.1 

17 

-10.5 

III 

nn 

--CDEFGH- 

2 

50 

4.9 

817 

923.0 

14.4 

16 

-11.0 

178 

3.0 

— CDEFGH- 

2 

52 

4.7 

823 

922.3 

14.9 

15 

-11.5 

III 

//// 

— CDEFGH- 

2 

54 

4.6 

827 

921.9 

15.0 

14 

-12.2 

III 

//// 

— CDEFGH- 

2 

56 

4.4 

833 

921.2 

15.1 

14 

-12.1 

III 

//// 

— CDEFGH- 

2 

58 

4.2 

841 

920.3 

15.1 

13 

-13.1 

III 

//// 

--CDEFGH- 

3 

0 

4.1 

843 

920.1 

15.2 

13 

-13.0 

181 

3.3 

--CDEFGH- 

3 

2 

3.9 

850 

919.4 

15.2 

13 

-13.0 

III 

llll 

--CDEFGH- 

3 

4 

4.0 

858 

918.5 

15.2 

12 

-14.0 

III 

llll 

— CDEFGH- 

3 

6 

3.9 

860 

918.3 

15.2 

12 

-14.0 

III 

llll 

--CDEFGH- 

3 

8 

3.5 

864 

917.8 

15.1 

12 

-14.1 

III 

llll 

— CDEFGH- 

3 

10 

3.5 

872 

916.9 

15.1 

12 

-14.1 

183 

3.7 

--CDEFGH- 

3 

12 

3.2 

880 

916.0 

15.1 

12 

-14.1 

III 

llll 

--CDEFGH- 

3 

14 

3.2 

885 

915.6 

15.1 

12 

-14.1 

III 

llll 

--CDEFGH- 

3 

16 

3.0 

889 

915.2 

15.1 

12 

-14.1 

III 

llll 

— CDEFGH- 

3 

18 

2.9 

897 

914.3 

15.1 

12 

-14.1 

III 

llll 

--CDEFGH- 

3 

20 

2.9 

903 

913.6 

15.2 

12 

-14.0 

186 

4.1 

— CDEFGH- 

3 

22 

2.9 

909 

912.9 

15.3 

12 

-13.9 

III 

llll 

--CDEFGH- 

3 

24 

2.9 

913 

912.5 

15.3 

12 

-13.9 

III 

llll 

— CDEFGH- 

3 

26 

2.7 

915 

912.3 

15.3 

12 

-13.9 

III 

llll 

--CDEFGH- 

3 

28 

2.9 

928 

910.9 

15.4 

12 

-13.8 

III 

llll 

--CDEFGH- 

3 

30 

3.0 

932 

910.5 

15.4 

11 

-14.9 

187 

4.6 

--CDEFGH- 

3 

32 

3.0 

940 

909.6 

15.5 

11 

-14.8 

III 

llll 

--CDEFGH- 

3 

34 

2.8 

942 

909.4 

15.6 

10 

-15.9 

III 

llll 

— CDEFGH- 

3 

36 

3.0 

951 

908.5 

15.7 

10 

-15.8 

III 

llll 

--CDEFGH- 

3 

38 

3.1 

957 

907.8 

15.8 

9 

-17.0 

III 

llll 

--CDEFGH- 

3 

40 

3.0 

961 

907.4 

15.9 

8 

-18.3 

189 

4.9 

— CDEFGH- 

3 

42 

3.0 

969 

906.5 

15.9 

8 

-18.3 

III 

llll 

— CDEFGH- 

3 

44 

3.0 

975 

905.8 

16.1 

7 

-19.7 

III 

llll 

--CDEFGH- 

3 

46 

3.1 

982 

905.2 

16.2 

7 

-19.6 

III 

llll 

--CDEFGH- 

3 

48 

3.1 

990 

904.3 

16.3 

7 

-19.6 

III 

llll 

--CDEFGH- 

3 

50 

3.0 

992 

904.1 

16.5 

7 

-19.4 

193 

5.0 

--CDEFGH- 

3 

52 

3.0 

998 

903.4 

16.5 

7 

-19.4 

III 

llll 

--CDEFGH- 

3 

54 

3.1 

1006 

902.5 

16.7 

7 

-19.3 

III 

llll 

--CDEFGH- 

3 

56 

3.3 

1013 

901.8 

16.7 

7 

-19.3 

III 

llll 

--CDEFGH- 

3 

58 

3.0 

1019 

901.2 

16.8 

7 

-19.2 

III 

llll 

--CDEFGH- 

4 

0 

3.2 

1027 

900.3 

16.8 

7 

-19.2 

196 

5.4 

--CDEFGH- 

4 

2 

3.0 

1031 

899.9 

16.9 

7 

-19.1 

III 

llll 

— CDEFGH- 

4 

4 

3.3 

1040 

899.0 

16.9 

7 

-19.1 

III 

llll 

— CDEFGH- 

4 

6 

3.0 

1042 

898.8 

16.9 

7 

-19.1 

III 

llll 

— CDEFGH- 

4 

8 

3.2 

1052 

897.7 

16.9 

7 

-19.1 

III 

llll 

--CDEFGH- 

4 

10 

3.2 

1056 

897.2 

16.9 

7 

-19.1 

197 

5.5 

--CDEFGH- 

4 

12 

3.0 

1058 

897.0 

17.0 

7 

-19.0 

III 

llll 

--CDEFGH- 

4 

14 

3.3 

1073 

895.5 

16.9 

7 

-19.1 

III 

llll 

--CDEFGH- 

4 

16 

3.0 

1073 

895.5 

16.8 

7 

-19.2 

III 

llll 

— CDEFGH- 

4 

18 

3.0 

1081 

894.6 

16.8 

7 

-19.2 

III 

llll 

--CDEFGH- 

4 

20 

3.1 

1085 

894.2 

16.8 

7 

-19.2 

198 

5.8 

— CDEFGH- 

4 

22 

3.1 

1091 

893.5 

16.8 

7 

-19.2 

III 

llll 

--CDEFGH- 

4 

24 

3.2 

1102 

892.4 

16.8 

7 

-19.2 

III 

llll 

— CDEFGH- 

4 

26 

3.1 

1106 

892.0 

16.8 

7 

-19.2 

III 

llll 

--CDEFGH- 

125 


126 


APPENDIX  F 


CONVERSION  OF  TEMPERATURE  TO  IN  BAND  RADIANCE 


%  Program  name :  rad . m 

%  Conversion  of  the  temperature  difference 
%  into  in-band  background  radiance 

clear; 

Tmin=238.;  %  Initial  temperature 


dT=  0.01;  %  Step 

n=2  00;  %  Number 

disp ( ' Temperature 

h=6 . 626e-34 ;  % 

c=2.998e8;  % 

k=l . 382e-23  ;  % 

a=8e-6; 

b=12 . 05e-6 ; 

dx=0 . 001e-6; 

for  T=Tmin : dT : Tmin+n 

sf=0 . 0 ; 

for  x=a:dx:b-dx 


of  values 

Radiance ' ) ; 

Planck's  constant 
Speed  of  light 
Boltzmann  constant 

dT 

%  Calculation  loop 


xm=x+dx/2 ; 

fl=  (2*h*c/N2)  /  (xm"5)  ; 
f2=l/ (exp ( (h*c) / (xm*k*T) ) -1) ; 
fx=f l*f 2 ; 
sf =sf +fx*dx; 

end; 

disp ( [T, sf ] ) ; 
end; 
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APPENDIX  G 


CALCULATION  OF  APPARENT  TEMPERATURE  DIFFERENCE  ATapp 

For  the  computation  of  the  effective  temperature  difference,  the  target  and 
background  temperatures  were  initially  converted  to  in-band  radiances  by  application  of 
the  Planck  function: 

N  -  £-^lM 

S  tr 

Tt 

where: 

£  =  emissivity 

a  =  Stefan-Boltzmann  constant  =  5.6697*  10'12  [W/m2K4] 

T  =  source  temperature  [K] 

Ag  =  fraction  of  the  thermally  emitted  radiance  lying  in  the  band  from  \  to  X2  at 
source  temperature  T 

So  it  is:  qi=/t,T,  qi=2jT  and  A^=q2-qi 

With  the  use  of  SEARAD  the  attenuated  radiance  and  path  radiance  were  calculated 
for  the  desired  range.  The  attenuation  from  the  polarizing  filter  was  also  calculated  and 
the  final  amount  of  energy  reaching  the  sensor  was  thus  obtained.  Then  by  using  the 
Planck  function,  the  radiance  was  converted  back  to  temperature  difference,  the  apparent 
temperature  difference  for  the  specific  range  and  geometry. 
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APPENDIX  H 


CALCULATION  OF  MRT  AND  MTF 

%  Program  mrt.m 
%  Calculation  of  MRT  and  MTF 
%  Lt  S.Lagaras  H.N. 

F0Vx=119 . 72959  ;  %  In  mRad  (=6.86  degrees) 

FOVy=90. 05899188;  %  In  mRad  (=5.16  degrees) 

SNRT=2.5;  %  Signal  to  Noise  Ratio 

%  NET=0.125;  %  Noise  Equivalent  Temperature  (if  you  want 

to  use  this  deactivate  other  NET) 

dx=0 . 1 ;  %  in  scan  detector  subtence  in  mrad 

dy=0.1;  %  cross  scan  detector  subtence  in  mrad 

L=7 ;  %  Length  to  width  ratio  of  the  bar  7 : 1 

Te=0.1;  %  Eye  integration  time  in  seconds 

Nss=l;  %  Serial  scan  ratio 

Nos=l ;  %  Overscan  ratio 

Fr=30;  %  Frame  rate  in  Hz  (1/sec) 

D=0 . 1 ;  %  Aperture  diameter  in  meters 

wl=10*10A (-6) ;  %  Approximation  of  the  radiation 

wavelength 

1=41.5;  %  Length  of  the  ship  in  meters 

w=9.75;  %  Width  of  the  ship  in  meters 

h=8.8;  %  Height  of  the  ship  in  meters 

%%%%%%  Calculation  of  back  end  resolution  (Shumaker  7-81) 
%rb=0.095;  %  Back  end  resolution  in  mRad 

td=6*10A ( -5 ) ;  %  Dwell  time  in  seconds 

f star=30*10A3 ;  %  Freq  at  which  the  detector  respone  is  down 
3db 

%  from  its  low  freq  response 

%  Detector  electronics  resolution  (rf) 

Vs=dx/td;  %  Angular  scan  velocity  in  mrad/ sec 
rf=Vs/ (fstar*pi) ; 

%  Preamplifier  resolution  (rp) 

DFp=8 . 5*10 A6 ;  %  Half  power  frequency  response  of  the 

preamplifier  in  HZ 
rp=0 . 32*pi*Vs/ (2*DFp) ; 
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%  Multiplexer  resolution  (rm) 

Ns=2;  %  #  of  multiplexer  samples  per  detector  dwell  time 

rm=dx/Ns ; 

%  Display  resolution  (rd) 

HD=0.2;  %  Display  height  in  meters 

dss=3 . 43*10A (-4) ;  %  Display  spot  size  in  meters 

rd=l . 9*dss*F0Vy/HD; 

%  Eye  resolution  (re) 

%  Assume  display  brightness  of  5  fl  so  that  re=1.2/M 
(Shumaker  pg  7-92) 

Dv=0.5;  %  Display  viewing  distance  in  meters 

M= (2*10A3*atan (HD/ (2*Dv) ) ) /FOVy;  %  System  magnification 
re=l . 2/M; 

%  Image  motion  resolution  (rs) 

AR=0.015;  %  Amplitude  of  residual  image  motion  in  mRad 

rs=3 . 55*AR;  %  (Shumaker  pg  7-94) 

rb=sqrt (rf A2+rpA2+rmA2+rdA2+reA2+rsA2 )  %  Back  end 

resolution  in  mrad 


%%%%%%  Variables  used  for  NET  calculation 


to=0 .70; 
Nd=180 ; 

Nsc=0 .15', 
Ds=4*10A10 ; 
DNDT=  6 . 2  7  * 1 0  A 
f  =  0 . 2  ; 


%  Transmission  of  the  optics 
%  Number  of  detectors 
%  Scan  efficiency 

%  Specific  detectivity  in  [cm  Hzl/2  Watt-1] 
(-5);  %  Derivative  of  Plank's  equation 

%  Focal  length  of  the  lens  in  meters 


N=13;  %  Number  of  bars  in  the  Bar  chart 

a=0  ; 

R=0. 0001:100:29000. 0001;  %  Range  in 

meters 

b=atan (1000 . /R) ; 

A= ( l*w*sin (b) ) + (h*w*cos (b) *sin(a) ) + (h*l*cos (b) *cos (a) ) ; 
At=abs (A) ; 

%%%%%%  At=214 . 4513 ;  This  is  the  value  used  in  other 

thesis 
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Dc=At.A(0.5)  ;  %  Critical  dimension 

X=Dc/N  ; 


fc=D/ (1000*wl) ;  %  Cut  off  frequency  in 

cycles /milliradians 

fs=R. / (2000*X) ;  %  Fundamental  spatial  freq  characterizing 

the  target  in  cycl/mrad 
m=f s . /fc ; 

Ho=  (2 /pi)  *  (acos  (m)  -  (m.  *  (1-m.  A2 )  .''0.5)); 

%Hd=sin (pi*dx*f s) / (pi*dx*fs) ; 

MTF=Ho; 

NET= ( 2  0  *  f  * (FOVx*FOVy*Fr*Nos*Nss ) *0 . 5) / (to* (pi*Nd*Nsc) A0 . 5*D 
A2*dx*dy*Ds*DNDT) ; 

Rx= ( 1+ (2*f s . *rb)  . A2 )  . A  (-0 . 5) ;  %  Shumaker  pg  8-60 

MRT=2  *SNRT*NET . * (Rx . A0 . 5 )  . *  (  ( (fs.  A2)  .*dx*dy/L)  . A0 .5)  .* ( (Te* 
Fr*Nos*Nss) A (-0 . 5) ) . /MTF ; 


figure (1) 
plot ( R , MRT ) ; 
axis( [0,29000,0,25] ) 
grid; 

ylabel (' Temperature  difference  in  Celcius ' ) ; 
xlabel ( ' Range  in  meters ' ) ; 

title ('MRT  vs  Range  for  Identification  N=13  (Johnson)'); 


figure (2 ) 

plot (m,MTF, 'r' ) ; 

grid; 

axis ([0,1. 2, 0,1]) 
ylabel ( 'MTF' ) ; 
xlabel ( ' f s/fc '  )  ; 
title ('MTF  vs  fs/fc'); 

figure (3 ) 

plot (fs,MTF, 'r ' ) ; 

grid; 

ylabel ( ' MTF ' ) ; 
xlabel ( ' f s ' )  ; 
title ('MTF  vs  f s ' )  ; 
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APPENDIX  I 
IDL  OUTPUTS 

A  number  of  IDL  programs  originally  created  by  M.C.  Pontes  and  later  modified  to 
fit  the  needs  of  this  work  were  used  to  obtain  the  temperature  information  of  R/V  POINT 
SUR  from  the  reordered  images.  The  programs  used  are  available  from  professor  A.W. 
Cooper.  In  the  following  image  is  presented  the  DDL  output  data  in  a  form  of  an  image 
that  was  used  for  this  work.  The  final  ship  temperature  was  finally  obtained  through  the 
average  temperature  of  the  two  thermistors  in  the  hot  spot  area. 


IDL  output  image. 
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APPENDIX  J 


EXCEL  SPREADSHEET  OUTPUTS 

The  following  table  is  part  of  the  EXCEL  spreadsheet  that  was  used  for  calculation  to 
the  apparent  temperature  difference  for  the  unpolarized  case  for  the  100m  scenario. 


File 

Range 

(KM) 

Zenith 

Angle 

Emissivity 

Hot  Spot 
Temp.(K) 

Ship  Radiance 
(0  Km) 

Transmittance 

Path  Radiance 

Ship+Path 

Radiance 

Ship+Path 

Temperature 

Background 

Radiance 

Background 

Temperature 

Apparent 

DT 

1 

0.107 

159.562 

0.9875 

300.23 

37.023798 

0.962 

1.03394 

36.64713177 

296.03 

32.47968 

288.97 

7.06 

15 

0.155 

130 

0.9748 

300.23 

37.023798 

0.952 

1.33149 

36.55963427 

295.88 

32.16871 

288.42 

7.46 

18 

0.199 

120 

0.9506 

300.23 

37.023798 

0.943 

1.58085 

36.48688722 

295.76 

31.63351 

21 

0.306 

109 

300.23 

37.023798 

0.924 

2.12944 

36.33572743 

295.52 

29.98601 

11.04 

24 

0.717 

98 

300.23 

37.023798 

0.864 

3.88881 

35.89218142 

279.05 

15.74 

95 

0.6856 

300.23 

5.43071 

35.53846294 

294.2 

27.00474 

278.78 

15.42 

92 

0.6236 

300.23 

10.28435 

34.54974753 

292.54 

27.34466 

279.45 

13.09 

91.5 

0.6147 

300.23 

37.023798 

0.586 

12.44314 

34.1427883 

291.85 

27.59885 

279.95 

11.9 

271 

5.307 

91.1 

0.6093 

300.23 

37.023798 

0.5 

15.13976 

33.65906401 

291.02 

27.96645 

10.36 

299 

6.551 

mm 

0.6075 

300.23 

37.023798 

0.439 

17.07819 

33.31682802 

■BEES 

28.26909 

9.19 

298 

7.43 

90.8 

0.6069 

300.23 

37.023798 

0.401 

18.2796 

33.1076313 

290.07 

28.44433 

281.58 

8.49 

297 

8.595 

0.6064 

300.23 

37.023798 

289.65 

28.68894 

282.04 

7.61 

296 

10.229 

0.6062 

300.23 

37.023798 

21.37795 

32.57394667 

289.14 

28.9886 

282.61 

6.53 

295 

12.723 

90.5 

0.6062 

HSUS] 

37.023798 

0.237 

23.43201 

32.21775738 

288.51 

29.36314 

283.32 

5.19 

2—13 

15.027 

90.44 

0.6063 

HESSES 

37.023798 

0.191 

24.9017 

31.95843599 

288.05 

29.63658 

283.83 

4.22 

2—11 

16.035 

90.42 

0.6063 

300.23 

37.023798 

6.173 

25.44248 

31.86240666 

287.88 

29.73818 

284.02 

3.86 

294 

17.225 

90.4 

0.6064 

300.23 

26.01343 

31.75582115 

287.69 

29.84631| 

284.22" 

3.47 

2-8 

18.668 

90.38 

0.6065 

300.23 

37.023798 

0.136 

287.5 

29.96226 

284.44 

3.06 
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APPENDIX  K 

SPECTRAL  RADIANCE  OF  A  BLACK  BODY  AT  DIFFERENT 

TEMPERATURES 

In  the  following  table  a  plot  of  the  spectral  radiances  for  a  black  body  at  303.3  K 
and  at  290  K  is  presented.  The  purpose  of  this  plot  is  to  evaluate  the  spectral 
characteristics  of  the  target  and  background  as  they  were  used  for  the  purposes  of  the 
thesis.  As  a  result  the  spectral  characteristics  of  the  target  and  background  do  not  seem  to 
have  any  effect  on  the  fact  that  the  Apparent  Temperature  Difference  for  the  polarized 
case  goes  negative  after  a  certain  range. 

This  plot  shows  the  distribution  of  the  radiance  with  wavelength  for  sub-bands 
across  8-12  fJm  band,  for  gray  body  sources  at  290  K  (sea  temperature)  and  303.3  K  (the 
ship  hot  spot  temperature).  The  wavelength  extend  of  the  sub-bands  is  not  constant.  From 
these  emitted  sub-band  radiances  the  available  radiances  at  25  km  were  computed  using 
the  sub-band  transmission  values  shown  in  Figure  31.  The  effective  band  transmittance  is 
then  available  as  the  ratio  of  the  total  radiance  at  25  km  to  total  emitted  radiance.  This 
quantity  computed  as  .0756  showed  no  difference  between  the  ship  and  background 
temperatures.  However,  it  differs  significantly  from  the  LOWTRAN  band  averaged 
transmittance  of  0.14.  It  should  be  noted  that  this  comparison  is  strictly  between  the  sea 
surface  black-body  emission  and  the  ship  hot  spot  temperature  black  body  emission. 
Reflected  sky  radiance  is  not  included  in  the  plot  numbers. 
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S  pectral  R adiance  at  T  =  303. 3| 
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